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Abstract 


The  influence  of  changes  in  discharge  and  season  on  the  concentration,  partitioning, 
and  export  of  a  limiting  nutrient,  phosphorus,  was  quantified  for  two  streams  in  central 
Alberta  from  April  to  October,  1983.  The  streams.  Two  Creek  and  Sakwatamau  River,  flow 
over  glacial  till  overlying  sedimentary  bedrock  and  are  located  on  coniferous  watersheds. 
Phosphorus  (P)  was  quantified  in  three  size  fractions:  1)  total  dissolved  phosphorus  (TDP; 
<0.45  wm),  2)  fine  particulate  phosphorus  (FPP;  0.45  urn  -  1  mm),  and  3)  coarse  particulate 
phosphorus  (CPP;  >1  mm).  Total  phosphorus  (TP)  is  the  sum  of  TDP  and  FPP. 

During  1983,  13.04  mg  P/m2  of  watershed  was  transported  past  the  Two  Creek 
sampling  station  and  7.49  mg  P/m2  was  transported  past  the  Sakwatamau  River  sampling 
station;  72%  of  the  P  load  was  in  the  FPP  fraction  for  both  streams.  Phosphorus  export  was 
strongly  pulsed  during  storm  events;  in  Two  Creek,  42%  of  TP  was  exported  with  19%  of  the 
water  in  2%  of  the  time.  Export  of  FPP  was  more  strongly  pulsed  than  was  export  of  TDP;  in 
Two  Creek  75%  of  FPP  was  transported  within  the  same  time  period  as  47%  of  TDP  transport. 
Surprisingly,  spring  runoff  was  not  the  most  significant  period  for  P  export.  In  fact,  65%  of 
annual  P  loading  was  transported  in  June  and  29%  in  July  in  Two  Creek.  This  was  probably 
due  to  low  winter  snowfall  relative  to  eastern  watersheds  and  heavy  rainstorms  with  overland 
flow  during  late  June  and  early  July.  Coarse  particulate  P  contributed  minimally  to  P  export 
(0.4  mg/m2  and  0.02  mg/m2  for  Two  Creek  and  Sakwatamau  River,  respectively). 

Models  to  predict  P  export  were  developed  for  the  study  streams.  In  these  streams, 
discharge  and  FPP  concentration  were  positively  correlated.  A  nonlinear  model  to  predict  FPP 
concentration  from  discharge,  which  included  rate  of  change  of  discharge,  was  significantly 
better  than  a  linear  model  for  the  study  streams.  There  were  no  significant  differences  in  the 
relationship  between  FPP  and  discharge  in  different  seasons.  Total  DP  was  also  positively 
correlated  to  discharge  in  both  streams  although  the  slope  of  TDP  concentration  on  discharge 
was  significantly  lower  than  the  slope  of  FPP  on  discharge. 
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Stream  sampling  methodologies  for  collecting  water  samples  were  evaluated.  When 
replicate  grab  samples  were  collected,  there  were  no  significant  differences  in  TP  and  TDP 
concentrations  in  either  stream.  Similarly,  when  samples  were  collected  across  the  width  of 
either  stream,  there  were  no  differences  in  FPP  or  TDP  concentrations  on  any  of  the  dates 
sampled.  Thus,  the  grab  sample  technique  used  in  this  study  was  adequate  to  estimate  P 
concentrations  in  streamwater.  Also,  there  were  no  diel  patterns  in  P  concentrations  in  Two 
Creek  and  Sakwatamau  River. 

The  results  of  longitudinal  transport  of  P  in  Two  Creek  suggest  seasonal  trends  in  size 
of  materials  in  transport  due  to  significant  differences  in  both  TDP  and  FPP  concentrations 
between  August  and  October.  However,  TDP  or  FPP  concentrations  did  not  change 
significantly  along  3390  m  of  Two  Creek  in  August  or  in  October  suggesting  no  longitudinal 
trends  in  size  of  materials  in  transport.  I  postulate  that  in  Two  Creek,  greater  transport  of 
particulate  material  during  summer  than  autumn  was  primarily  due  to  higher  flows  and  not 
invertebrate  processing. 

This  study  concurs  with  earlier  predictions  that  abiotic  factors,  such  as  watershed  land 
use,  bedrock  geology,  and  climatic  regimes,  all  influence  P  export  patterns.  However,  this  study 
suggests  that  P  export  coefficients  based  on  watershed  area  may  not  be  an  appropriate  export 
expression  for  forested  streams  overlying  sedimentary  bedrock. 
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I.  General  Introduction 


Transfers  of  matter  or  energy  involving  one  or  more  ecosystem  component  have  been 
studied  in  attempts  to  understand  stream  ecosystem  processes  and  functioning  (e.g.,  Cummins 
1974,  Likens  and  Bormann  1974,  MacMahon  et  al.  1981).  Phosphorus  (e.g.,  Hobbie  and  Likens 
1973,  Cahill  et  al.  1974,  Johnson  et  al.  1976,  Meyer  and  Likens  1979,  Rigler  1979,  Hill  1982) 
and  carbon  (e.g.,  Hobbie  and  Likens  1973,  Wetzel  and  Manny  1977,  Tate  and  Meyer  1981, 
Wallace  et  al.  1982)  are  two  frequently  studied  components.  Phosphorus  (P)  is  relatively 
insoluble  and  unlike  other  nutrients  such  as  carbon  and  nitrogen,  P  has  no  gaseous  form. 
Although  P  is  found  in  low  levels  in  natural  water,  it  is  essential  for  algal  and  higher  plant 
growth.  The  addition  of  P  to  surface  waters  can  enhance  autotrophic  production  and  rates  of 
detrital  processing  (Schindler  1977,  Stockner  and  Shortreed  1978,  Elwood  et  al.  1981).  For 
these  reasons,  P  is  often  considered  to  limit  production  in  undisturbed  aquatic  systems.  Thus, 
by  defining  the  energy  pathways  of  a  fundamental  component  such  as  P,  key  processes  in  a 
stream  ecosystem  may  be  predicted  from  deterministic  and  stochastic  changes  in  the  ecosystem. 

To  understand  the  energy  dynamics  of  an  essential  component,  it  is  necessary  to  define 
the  inputs  and  outputs  to  the  ecosystem.  Meyer  and  Likens  (1979)  present  the  most  rigorous 
example  of  this  approach  at  Bear  Brook,  an  undisturbed  headwater  stream  in  the  Hubbard 
Brook  Experimental  Forest,  New  Hampshire.  A  P  mass  balance  was  constructed  for  Bear 
Brook  for  a  13 -yr  period  based  on  an  empirical  model  of  the  annual  budget.  Nutrient  inputs 
were  quantified  in  three  size  fractions:  1)  total  dissolved  P  (<0.45  um),  2)  fine  particulate  P 
(0.45  um  -  1  mm),  and  3)  coarse  particulate  P  (>1  mm).  Phosphorus  inputs  to  Bear  Brook 
for  all  three  size  fractions  were  from  precipitation,  subsurface  waters,  tributaries,  and  leaf  litter 
falling  directly  into  the  stream.  The  important  exports  were  all  fluvial.  In  Bear  Brook,  P 
accumulated  and  was  processed  in  the  stream  ecosystem  on  most  days  of  the  year.  Export  was 
confined  to  short  periods  of  high  stream  discharge,  and  the  ratio  of  P  outputs  to  inputs  was 
directly  related  to  annual  streamflow.  The  net  effect  of  ecosystem  processes  at  Hubbard  Brook 
was  the  conversion  of  dissolved  P  and  coarse  particulate  P  into  fine  particulate  P,  which  was 
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then  transported  downstream. 

However,  most  studies  have  more  limited  objectives;  they  quantify  exports  to 
downstream  water  bodies  (e.g.,  Hobbie  and  Likens  1973,  Johnson  et  al.  1976,  Qngley  1976, 
Wartiovaara  1978,  Rigler  1979,  Hill  1981).  The  "export  coefficient"  approach  (mass  per  unit 
area  of  watershed)  is  used  to  study  the  effect  of  one  ecosystem  component  on  another,  such  as 
the  effect  of  land  use,  pollution,  and  nutrient  fertilization  on  nutrient  levels  in  streams  or  lakes 
(e.g.,  Ryden  et  al.  1973,  Dillon  and  Kirchner  1975,  Qmernik  1977,  Rast  and  Lee  1978,  Stockner 
and  Shortreed  1978),  and  has  influenced  the  development  of  nutrient  loading  calculations  for 
eutrophication  management.  This  approach  has  also  been  used  to  determine  the  effect  of 
clear-cutting  on  P  export  (e.g.,  Bormann  et  al.  1968,  Hornbeck  1975)  and  to  study  seasonal 
changes  in  loading.  Temporal  or  season -related  changes  in  the  ecosystem  (e.g.  snowmelt  and 
storm  events)  influence  P  export  and  can  vary  from  year-to-year  (Meyer  and  Likens  1979). 
Seasonal  variation  in  P  export  has  received  increased  attention  in  recent  years  because  natural 
ecosystem  variation  must  be  understood  before  variation  can  be  attributed  to  perturbations 
such  as  logging.  I  will  discuss  the  effects  of  temporal  changes  on  P  export  in  more  detail  later. 

Approaches  that  focus  on  nutrient  inputs  and  outputs  to  streams  give  information  only 
on  net  transfers  and  not  on  transformations  within  the  system.  For  example,  P  in  a 
decomposing  leaf  may  be  assimilated  into  fungal  hyphae;  both  compounds  are  a  similar  size 
fraction  (>1  mm)  even  though  they  are  very  different  compounds  (Meyer  1978).  While  Meyer 
(1978)  and  earlier  stream  biologists  largely  ignored  this  type  of  processing  in  their  ecosystem 
studies,  emphasis  has  recently  been  placed  on  understanding  physical,  chemical,  and  biological 
processes  that  control  material  cycling  (e.g.,  Newbold  et  al.  1982,  Minshall  et  al.  1983).  The 
cycling  concept  (termed  spiralling)  theorizes  that  the  productivity  and  persistence  of  stream 
ecosystems  depend  on  the  retention  and  recycling  of  nutient  and  energy  resources  from  the 
surrounding  watershed  as  they  are  transported  downstream  (Newbold  et  al.  1982).  Also,  the 
degree  of  retention  and  reutilization  of  nutrients  in  a  stream  is  associated  with  the  "tightness" 
of  the  spirals  (Webster  1975). 
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Newbold  and  his  colleagues  (1982)  at  the  Oak  Ridge  National  Laboratory  in  Tennessee 
constructed  a  steady -state  model  to  predict  spiralling  length  of  a  limiting  nutrient.  A  complete 
cycle  was  considered  to  consist  of  biotic  uptake  of  a  nutrient  atom  from  an  available  dissolved 
state,  subsequent  passage  through  the  food  chain,  and  ultimate  return  to  the  water  in  an 
available  dissolved  form.  Soluble  reactive  P  (SRP)  was  measured  as  an  estimate  of  biologically 
available  dissolved  P.  The  model  was  based  on  five  assumptions:  (1)  the  system  is  in 
steady -state,  (2)  microbial  growth  and  algal  photosynthesis  are  limited  through  a  concentration 
limitation  on  nutrient  uptake  kinetics,  (3)  all  available  P  is  within  the  SRP  form,  (4)  all  SRP  is 
available,  and  (5)  soluble  unreactive  P  is  not  important.  The  model  predicts  that  most 
downstream  transport  of  the  nutrient  occurs  in  an  unavailable  particulate  form  when  nutrient 
limitation  is  severe  but  when  "nutrient  limitation  is  moderate",  transport  in  the  dissolved  phase 
dominates  and  how  well  particles  are  transported  has  little  influence  on  spiralling  length. 
However,  the  terms  "moderate"  and  "severe"  nutrient  limitation  have  never  been  operationally 
defined.  The  nutrient  spiralling  model  is  important  theoretically  because  it  can  incorporate  the 
potential  role  of  invertebrate  consumers  on  nutrient  dynamics  and  could  be  elaborated  to 
incorporate  other  levels  of  the  food  web  (e.g.  fish).  In  a  recent  study,  radioactive  P  was 
released  into  a  woodland  stream  in  Tennessee  and  its  uptake  from  water  to  coarse  particulate 
organic  matter  (CPOM),  fine  particulate  organic  matter  (FPOM),  aufwuchs,  grazers,  etc.,  was 
measured  (Newbold  et  al.  1983).  In  this  study,  the  average  velocity  of  P  travelling  downstream 
was  10.4  m/d  and  P  completed  the  cycle  once  every  18.4  d.  The  average  downstream  distance 
associated  with  one  cycle,  defined  as  spiralling  length,  was  therefore  190  m  (10.4  m/d  X  18.4 
d).  Coarse  particulate  organic  matter  accounted  for  60%  of  the  uptake  although  FPOM 
transport  accounted  for  99%  of  the  particulate  turnover  length  due  to  a  slow  turnover  time  (99 
d).  As  well,  P  retention  within  the  stream  community  was  high;  Newbold  and  his  colleagues 
attribute  this  to  low  drift  rates  of  consumers  and  a  long  turnover  time  for  P  (152  d)  within  the 


consumer  community. 


- 


4 


Another  approach  for  defining  stream  ecosystem  processes  is  the  river  continuum 
concept.  It  was  developed  to  describe  gradual  changes  in  ecosystem  structure,  function,  and 
stability  along  the  length  of  a  stream  (Vannote  et  al.  1980).  This  concept  predicts  that  as 
stream  size  increases,  the  reduced  importance  of  allochthonous  organic  input  coincides  with  the 
enhanced  significance  of  autochthonous  primary  production  and  organic  transport  from 
upstream;  the  latter  is  reflected  by  an  increase  in  the  ratio  of  gross  primary  productivity  to 
community  respiration  (see  Vannote  et  al.  1980).  Unlike  the  nutrient  spiralling  model,  the  river 
continuum  concept  is  based  on  a  dynamic  equilibrium  whereby  producer  and  consumer 
communities  characteristic  of  a  river  reach  become  established  in  harmony  with  dynamic 
physical  conditions  of  the  channel.  However,  the  river  continuum  concept  is  too  general,  and 
deviations  from  expected  patterns  are  too  easily  explained  (e.g.  Minshall  et  al.  1983). 

It  is  apparent  that  stream  ecosystem  studies  are  elaborate  and  time  consuming.  Thus,  it 
is  important  to  have  well  defined  goals  within  budget  and  time  constraints.  A  decade  ago,  many 
studies  were  conducted  to  evaluate  the  effect  of  land  use  on  P  loading  in  different  biomes  (e.g. 
Dillon  and  Kirchner  1975,  Omernik  1977,  Rast  and  Lee  1978).  Most  of  these  projects  were 
government  directed  with  eutrophication -control  goals.  In  these  projects,  many  different 
streams  were  sampled  once  every  two  to  three  weeks.  Dillon  and  Kirchner  (1975)  listed  annual 
P  export  values  for  streams  situated  in  watersheds  with  different  geology  and  land  use  from  all 
over  the  world.  The  most  widely  used  theoretical  export  coefficients  for  nutrient  budgets  were 
based  on  data  from  34  watersheds  in  southern  Ontario  (Dillon  and  Kirchner  1975)  or  on  data 
from  all  across  the  United  States  (Rast  and  Lee  1978).  More  recently,  Meyer  and  Likens  (1979) 
and  Rigler  (1979)  demonstrated  that  most  P  export  values  reported  in  the  literature  are 
underestimates  of  actual  loading  because  the  data  are  derived  from  inadequate  sampling  regimes 
(2-3  wk  intervals)  which  neglected  short-term  changes  resulting  from  storm  runoff  (see  Table 
1.1).  Stevens  and  Smith  (1978)  made  a  comparison  between  loading  values  obtained  from  2-8  h 
sampling  intervals  versus  8  d  sampling  intervals  over  a  102  d  period.  They  found  particulate  P 
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was  underestimated  by  43%  with  8  d  as  compared  to  2-8  h  sampling  intervals;  elevated 
particulate  P  concentrations  during  short  periods  of  high  discharge  were  missed  with  the  8  d 
routine.  Also,  soluble  reactive  P  (potentially,  the  fraction  of  dissolved  phosphorus  that  is 
available  for  uptake)  was  overestimated  by  12%  with  the  8  d  routine  because  soluble  reactive  P 
concentration  decreased  at  high  discharge  due  to  dilution  by  rainwater.  A  similar  comparison 
was  done  for  a  moorland  stream  in  England;  annual  export  was  recalculated  as  if  the  samples 
had  been  only  collected  on  a  weekly  or  monthly  basis  (Rigler  1979).  Both  weekly  and  monthly 
export  calculations  were  significantly  less  than  'true  annual  export'  based  on  an  intensive 
sampling  program.  Thus,  loading  values  in  empirical  models  based  on  a  1-2  wk  sampling 
interval  are  poor  estimates  and  require  refinements. 

Regional  variations  in  general  ecosystem  characteristics  (e.g.  geology,  slope),  timing  of 
seasonal  events,  and  stochasticity  of  storms  result  in  a  variety  of  loading  patterns.  A  stream  in 
rural  New  York  (43%  agricultural  land)  exported  75%  of  annual  P  during  highest  flows  which 
covered  only  10%  of  the  time  (Johnson  et  al.  1976).  During  snowmelt  runoff,  50%  of  total 
annual  P  export  was  transported  from  two  contour -farmed  Iowa  watersheds  (Alberts  et  al. 
1978).  In  only  10  days,  67%  of  the  annual  P  export  was  transported  from  a  forested  stream  in 
New  Hampshire  (Meyer  and  Likens  1979).  A  large  portion  of  annual  stream  P  inputs  can  occur 
within  the  first  few  days  or  hours  of  spring  runoff  in  Arctic  streams  (de  March  1975,  H.E. 
Welch,  Freshwater  Institute,  pers.  comm.).  Thus,  P  loading  data  exhibit  strong  temporal  and 
seasonal  trends,  generally  with  peaks  during  snowmelt  and  storm  events. 

Regional  variations  exist  in  the  relationship  between  P  concentration  and  discharge. 
Increasing  discharge  resulted  in  marked  increases  in  total  P  concentration,  especially  in  the  fine 
particulate  fraction,  for:  (1)  a  forested  stream  (Bear  Brook)  on  the  Precambrian  Shield 
(Meyer  and  Likens  1979),  (2)  grazed  moorland  streams  of  Dartmoor,  England  (Rigler  1979), 
and  (3)  arctic  streams  entering  Char  Lake  (de  March  1975).  In  Bear  Brook  and  the  Dartmoor 
streams,  increases  in  P  concentration  during  floods  were  larger  during  ascending  than 
descending  stage  of  the  hydrograph,  and  P  increases  were  greater  in  summer  than  winter 
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(Meyer  and  Likens  1979,  Rigler  1979).  In  these  two  streams,  most  of  the  P  increase  during 
flooding  was  in  the  particulate  fraction;  Rigler  (1979)  postulated  that  these  increases  could  be 
explained  if  P  accumulated  on  the  stream  bed  during  periods  of  base  flow  and  was  resuspended 
as  particulate  P  during  floods.  In  Bear  Brook,  the  relationship  between  fine  particulate  P  and 
discharge  varied  with  season,  position  of  hydrograph,  and  watershed  area.  In  the  moorland 
streams,  dissolved  P  concentrations  increased  slightly  with  an  increase  in  discharge  whereas  in 
Bear  Brook,  there  was  no  significant  change  in  dissolved  P  concentration  when  discharge 
increased.  In  streams  draining  Ohio  woodlands  (Taylor  et  al.  1971)  and  in  two  rivers  draining 
disturbed  watersheds  in  southern  Alberta  (Bow  River  and  Oldman  River),  dissolved  P 
concentration  and  discharge  were  not  related  (Qngley  and  Blachford  1982). 

Studies  conducted  in  the  larger  rivers  in  western  Canada  have  concentrated  primarily  on 
longitudinal  trends  in  P  concentration;  year-to-year  variation  was  monitored,  and  effects  of 
dumping  urban  effluents  into  the  river  evaluated  (Paterson  and  Nursall  1975,  Tones  et  al.  1980, 
Ongley  and  Blachford  1982).  Studies  on  small  streams  entering  lakes  have  been  conducted  for 
Lake  Wabamun  (Mitchell  et  al.  1981),  Baptiste  Lake  (Trew  et  al.  1978),  Cold  Lake  (Trew  et 
al.  1981),  and  Nakamun  Lake  (Riley  and  Prepas  1984)  as  part  of  lake  nutrient  budgets.  Most 
streams  were  intermittent.  Trew  and  his  colleagues  at  Alberta  Environment  (1978)  concluded 
that  Rast  and  Lee's  (1978)  export  coefficients  were  adequate  for  streams  draining  into  Baptiste 
Lake.  However,  P  loading  data  for  both  Baptiste  Lake  and  Rast  and  Lee's  loading  coefficients 
were  based  on  inadequate  sampling  regimes. 

Many  gaps  exist  in  our  knowledge  of  patterns  in  stream  phosphorus.  Good  export 
coefficients  have  yet  to  be  developed  for  watersheds  on  glacial  till.  Loading  values  have  never 
been  presented  on  a  seasonal  basis  for  non -shield  streams  and  this  may  be  important  in  the 
design  of  monitoring  programs.  As  well,  the  relationship  between  P  concentration  and 
discharge  has  never  been  adequately  tested  for  undisturbed  streams  on  glacial  till.  The 
contribution  of  coarse  particulate  P  to  P  loading  has  been  investigated  for  only  one  stream. 
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Bear  Brook,  which  is  in  a  deciduous  forest  and  flows  over  igneous  bedrock  (Meyer  and  Likens 
1979).  Also,  the  standard  approach  for  sampling  stream  water,  the  grab  sampling  technique, 
has  never  been  tested  for  adequacy  in  estimating  P  concentration  in  streamwater.  In  other 
words,  horizontal  heterogeneity  in  stream  P  concentration  at  one  location  has  not  been 
evaluated.  In  this  study,  I  investigated  these  questions  in  two  undisturbed  streams  located  in  a 
boreal  mixed -wood  biome  in  central  Alberta;  I:  (1)  investigated  the  influence  of  season  on  P 
loading,  (2)  examined  the  influence  of  deterministic  (e.g.  spring  runoff)  and  stochastic  (e.g. 
storm  events)  changes  in  the  ecosystem  on  P  concentration  and  partitioning,  (3)  evaluated  the 
grab  sample  technique,  and  (4)  developed  a  model  to  predict  P  loading  in  headwater  streams 
flowing  over  3-5  m  of  glacial  till  overlying  sedimentary  bedrock. 
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TABLE  1.1.  Annual  phosphorus  (P)  loading  for  various  watersheds  in  north  temperate  zones.  Data  are  divided  into  two  groups  based  on  whether  the 
sampling  regime  was  sufficiently  frequent  (adequate)  or  not  sufficiently  frequent  (inadequate)  to  produce  accurate  loading  coefficients.  Location, 
geology,  and  land  use  are  indicated. 


12 


O  CO  ri  '"'I  0  0 

_!  rX  ^  r— t  ^ 


r- 


T“— C 

^  qn 

rs  T-H 


r-* 


"O 
T3  0> 
CD  JD 


0) 


°  iS 

U-H 


a> 

T3 

T3 

Vh 

CD 

CD 

3 

S3 

03 

Wh 

S3 

03 

Vm 

4— » 

3 

00 

00 

o 

N 

S's* 

4-. 

T3 

•d  T?  'd  T3 

T3  OQ 

<D 

<D 

<D 

CD 

CD 

0>  03 

— -> 

4— > 

■w> 

4— > 

•*->  \ 

C/3 

C/3 

(30 

C/3 

C/3  * 

0> 

CD 

CD 

<D 

0)  03 

V-. 

t-4 

W-. 

Wh 

W-.  t-4 

to 

O 

o 

o 

O 

O  =3 

Uh 

<4_ 

U-, 

<— 4  t-4 

*3 

<D 

S3 

03 

u. 

oo 


T3 

<D 

S3 

OJ 

t- 


01) 


>> 


C/3 

C/3 

03 

4— > 

G 

C/3 

C/3 

c/3 

=3 

=3 

CD 

*"0 

G 

G 

O 

O 

E 

D 

O 

O 

O 

CD 

CD 

CD 

CD 

<D 

c 

c 

*■5 

t-4 

G 

g 

c 

oo 

OO 

CD 

03 

OO 

oo 

oo 

oo 

>• 

on 

c 

K 

H 

D 

< 


Meyer  and  Likens  (1979)  Hubbard  Brook,  N.H.  igneous  forested 


II.  Seasonal  dynamics  of  phosphorus  partitioning  and  export  in  two  headwater  streams  in  central 

Alberta 


A.  Abstract 

In  1983,  phosphorus  (P)  export  was  quantified  for  two  undisturbed  headwater  streams 
in  central  Alberta,  Canada.  These  streams  flow  over  3-5  m  of  glacial  till  overlying  sedimentary 
bedrock  and  are  located  in  boreal  mixed -wood  watersheds.  In  these  streams,  the  influence  of 
changes  in  discharge  on  P  concentration  and  partitioning  was  examined  on  an  annual  and 
seasonal  basis  and  compared  with  streams  on  granitic  bedrock  with  disturbed  or  hardwood 
forested  watersheds.  The  relationship  between  P  concentration  and  discharge  was  used  to 
develop  models  to  predict  P  concentration  in  the  forested  Albertan  streams.  Unlike  previous 
studies  where  P  export  was  maximum  in  the  spring  or  autumn,  P  export  in  the  study  streams 
peaked  during  summer  storm  events;  68%  of  total  P  loading  was  transported  during  12  days  in 
early  summer.  In  the  study  streams,  hysteresis  effects  were  not  observed  during  the  early 
summer  storms  but  were  observed  during  later  summer  storms.  Phosphorus  export  was 
primarily  in  the  fine  particulate  fraction  (0.45  urn  -  1  mm;  72%).  Coarse  particulate  P  was  the 
least  important  fraction  (>  1  mm;  3%  and  0.3%  for  the  two  streams)  possibly  due  to  the 
absence  of  an  autumn  deciduous  leaf  fall.  There  were  strong  positive  correlations  between 
discharge  and  both  fine  particulate  P  and  total  dissolved  P  (<  0.45  um);  the  correlations  were 
nonlinear  and  linear,  respectively.  Models  presented  are  the  first  based  on  detailed  data  from 
watersheds  with  the  following  characteristics:  (1)  coniferous  trees,  (2)  3-5  m  of  glacial  till 
overlying  sedimentary  bedrock,  and  (3)  climatic  regime  where  summer  storms  are  a  dominant 
influence  on  P  export. 
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B.  Introduction 

A  decade  ago,  many  studies  were  conducted  to  evaluate  the  effect  of  land  use  on 
phosphorus  (P)  loading  in  different  biomes  (e.g.  Dillon  and  Kirchner  1975,  Rast  and  Lee 

1978) .  Most  of  these  projects  were  government  directed  with  eutrophication -control  goals.  In 
these  projects,  many  different  streams  were  sampled  once  every  two  to  three  weeks.  More 
recently,  Stevens  and  Smith  (1978),  Meyer  and  Likens  (1979),  and  Rigler  (1979)  demonstrated 
that  most  P  export  values  reported  in  the  literature  are  underestimates  of  actual  loading;  data 
derived  from  sampling  regimes  of  2 -to -3  wk  intervals  are  inadequate  because  they  miss 
dramatic  short-term  changes  during  storm  runoff.  For  example,  in  only  10  days,  67%  of  the 
annual  P  export  was  transported  from  a  forested  stream  in  New  Hampshire  (Meyer  and  Likens 

1979) . 

The  transport  of  P  within  different  size  fractions  and  the  relationship  with  discharge 
has  been  studied  for  two  groups  of  streams  flowing  over  granite:  the  headwater  streams  located 
in  deciduous  forest  at  Hubbard  Brook  (Meyer  and  Likens  1979)  and  the  grazed  moorland 
streams  in  Dartmoor,  England  (Rigler  1979).  In  the  Hubbard  Brook  and  Dartmoor  streams, 
increases  in  P  concentration  during  floods  were  larger  during  ascending  than  descending  stage 
of  the  hydrograph,  and  P  increases  were  greater  in  summer  than  in  winter.  Most  of  the  P 
increase  during  flooding  was  in  the  particulate  fraction;  Rigler  postulated  that  these  increases 
could  be  explained  if  P  accumulated  on  the  streambed  during  periods  of  base  flow  and  was 
resuspended  as  particulate  P  during  floods.  In  Hubbard  Brook,  the  relationship  between  fine 
particulate  P  (FPP;  0.45  um  -  1  mm)  and  discharge  varied  with  season,  position  of 
hydrograph,  and  watershed  area.  In  the  moorland  streams,  total  dissolved  P  (TDP;  <  0.45 
um)  concentrations  increased  slightly  with  an  increase  in  discharge  whereas  in  Hubbard  Brook, 
there  was  no  significant  change  in  TDP  when  discharge  increased.  As  well,  the  contribution  of 
coarse  particulate  P  (CPP;  >  1  mm)  to  P  loading  and  the  variation  in  P  concentrations 
between  summer  and  winter  have  been  investigated  for  only  one  stream  (Meyer  and  Likens 
1979).  Very  little  is  known  about  what  to  expect  in  different  stream  types  such  as  those 
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underlain  by  glacial  till  rather  than  granite  bedrock,  where  groundwater  seepage  may  be 
important,  where  seasonal  discharge  patterns  vary  from  the  classical  spring  and  autumn  peaks, 
and  where  vegetation  is  not  characterized  by  deciduous  trees  with  a  potentially  heavy  P  input  at 
autumn  leaf  fall. 

The  standard  approach  for  sampling  stream  water,  the  grab  sample  technique,  has  never 
been  tested  for  adequacy  in  estimating  P  concentration  in  stream  water.  In  other  words, 
horizontal  heterogeneity  in  stream  P  concentration  at  one  location  has  not  been  adequately 
evaluated.  As  well,  diel  patterning  may  influence  P  loading.  Although  Meyer  and  Likens  (1979) 
reported  no  diel  patterns  in  TDP  concentration,  de  March  (1975)  reported  diurnal  changes  in 
flow,  TDP,  and  FPP  concentrations  for  streams  entering  Char  Lake.  If  diel  patterns  exist, 
sampling  regimes  should  incorporate  samples  taken  at  different  times  in  the  24 -h  cycle. 

There  are  no  published  data  on  seasonal  variation  in  P  loading  and  partitioning  in 
streams  located  on  glacial  till  in  watersheds  with  primarily  coniferous  trees.  This  information  is 
necessary  both  to  develop  sampling  routines  to  monitor  nutrient  export  from  similar  streams  in 
boreal  forested  regions  and  to  develop  a  broader  perspective  on  seasonal  patterns  in  nutrient 
loading  and  P  partitioning  in  streams.  Thus,  I  investigated  the  following  questions  in  two 
undisturbed,  headwater  streams  located  on  glacial  till  in  a  boreal  mixed -wood  biome  in  central 
Alberta:  (1)  the  influence  of  season  on  P  export,  (2)  the  influence  of  deterministic  (e.g.  spring 
runoff)  and  stochastic  (e.g.  storm  events)  changes  in  the  ecosystem  on  P  concentration  and 
partitioning,  and  (3)  the  influence  of  the  grab  sample  technique  and  possible  diel  patterns  on 
sampling  programs.  I  also  developed  models  to  predict  P  loading  in  these  forested  headwater 


streams. 
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C.  Materials  and  Methods 

Field  Sampling-  Field  work  was  conducted  from  stream  ice-off  (April  22)  to  freeze-up 
(October  25)  during  1983.  Two  streams.  Two  Creek  and  Sakwatamau  River,  were  sampled 
twice  daily  and  more  frequently  during  storm  events  except  for  sixteen  2-d  periods  when 
samples  were  being  analyzed  in  the  laboratory.  Air  and  water  temperatures  were  measured  with 
a  hand  thermometer.  A  Price  A  A  current  meter  was  used  to  measure  velocity  and  water  depth 
at  1.0  or  0.5  m  intervals  across  the  stream  width.  Water  level  was  recorded  with  continuous 
recorders  as  well  as  from  a  depth  gauge  placed  in  each  stream  in  April.  On  each  visit,  a  grab 
sample  of  stream  water  was  collected  from  the  middle  of  each  stream  in  1-L  Nalgene  bottles  for 
P  analysis.  Occasionally,  grab  samples  of  stream  water  were  collected  for  analysis  of  other 
physical  and  chemical  parameters  such  as  turbidity,  colour,  and  conductivity.  To  estimate 
transport  of  CPP  (>1  mm),  coarse  particulate  matter  (CPM)  was  collected  from  a  1-mm 
mesh  net  (35  X  40  cm)  placed  in  each  stream  such  that  the  entire  water  column  was  sampled. 
To  quantify  the  volume  of  water  passing  through  the  nets,  water  velocity  was  measured  1  m  in 
front  of  each  net  at  the  time  of  sampling.  Material  was  collected  daily  from  the  nets  from  June 
to  October  and  later  analyzed  for  P  content. 

To  test  if  one  grab  sample  was  adequate  to  quantify  P  transported  downstream, 
replicate  grab  samples  were  taken  monthly  from  both  streams.  To  test  for  spatial  heterogeneity 
across  the  stream  width,  grab  samples  were  taken  at  2-m  intervals  across  the  width  of  the 
stream  and  analyzed  for  total  phosphorus  (TP)  and  TDP  concentrations.  To  test  for  diel 
variation  in  FPP  and  TDP  concentrations,  grab  samples  of  stream  water  were  collected  every  4 
h  for  24  h.  Samples  for  diel  analyses  were  collected  on  five  dates  for  Two  Creek  and  four  dates 
for  Sakwatamau  River  throughout  the  sampling  year. 

Daily  precipitation  data  were  obtained  from  an  Alberta  Forest  Service  fire  lookout 
station  located  on  the  northern  border  of  Two  Creek  watershed.  Discharge  data  for 
Sakwatamau  River  were  compared  to  discharge  data  collected  by  Environment  Canada  since 
1973  at  the  mouth  of  Sakwatamau  River  (Environment  Canada  1973-1983). 
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Heavy  rainstorms  flooded  the  road  to  Sakwatamau  River  from  July  4-19  and  the 
stream  could  not  be  sampled.  Consequently,  Sakwatamau  River  loading  data  were  interpolated 
from  a  regression  of  Sakwatamau  River  P  loading  data  on  Two  Creek  P  loading  data  for  the 
summer  months  excluding  this  16 -d  period  (r=0.92).  As  well,  Sakwatamau  River  data  were 
not  used  to  analyze  the  influence  of  July  storms  on  P  concentration  and  partitioning.  During 
the  heavy  rain  storms  in  July,  water  depth  increased  in  both  streams  such  that  normal  sampling 
procedures  were  disrupted;  flow  velocity  was  calculated  by  multiplying  surface  velocity  by  a 
factor  of  0.8  (Gray  1970).  To  collect  CPM,  1-mm  mesh  nets  were  held  from  the  bank  in  the 
water  for  3  to  5  min. 

The  two  study  streams  were  frozen  to  the  bottom  from  November  to  April  and 
consequently,  there  was  no  export  during  this  period. 

Laboratory  Methods-  Total  P  was  determined  in  triplicate  (standard  error  always  <  7% 
of  P  value)  on  50 -mL  subsamples  filtered  through  a  1-mm  mesh  net  and  analyzed  by  the 
potassium  persulfate  technique  (Prepas  and  Rigler  1982).  Water  to  be  analyzed  for  TDP  was 
filtered  through  a  prerinsed  0.45-um  Millipore  HAWP  membrane  filter  and  analyzed  as  for  TP. 
Fine  particulate  P  was  calculated  as  the  difference  between  TP  and  TDP.  To  determine  the  size 
of  materials  within  the  fine  particulate  size  fraction,  FPP  was  subdivided  further  (1  mm -80 
um,  80-64  urn,  64-10  um,  10-5  um,  5-0.45  wm)  on  eight  dates  from  August  to  October  as  well 
as  on  three  dates  in  the  spring  of  1984.  Coarse  particulate  matter  was  dried  for  24  h  at  105°C 
and  weighed.  Approximately  monthly,  six  1-g  subsamples  of  dried  CPM  were  analyzed  for  P 
content  and  these  data  were  used  to  convert  daily  CPM  dry  weights  to  daily  CPP  values. 

Total  Kjeldahl  nitrogen  was  analyzed  by  the  sulfuric  acid -copper  sulfate  method  of 
D’Elia  et  al.  (1977)  as  modified  by  Prepas  and  Trew  (1983).  Sulfate  analyses  were  done  on 
samples  filtered  through  a  prewashed  Whatman  GF/C  filter  by  the  gravimetric  method  with 
ignition  of  the  residue  (American  Public  Health  Association  1975).  Chloride  samples  were 
pretreated  as  sulfate  samples  and  analyzed  by  the  mercuric -nitrate  method  (American  Public 
Health  Association  1975).  Sodium  and  potassium  were  measured  with  flame  photometry 
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methods  (American  Public  Health  Association  1975).  Calcium  and  magnesium  were  analyzed 
with  flame  emission  on  a  Jarrel  Ash  (Model  JA82-270)  atomic  absorption  spectrophotometer 
(American  Public  Health  Association  1975)  except  no  lanthanum  was  added.  Alkalinity  was 
determined  by  the  potentiometric  titration  method  of  Environment  Canada  (1979).  A  Fisher 
Accumet  Digital  pH  meter  model  520  was  used  to  measure  pH  and  to  determine  the  end  points 
of  the  alkalinity  titration.  Colour  was  measured  visually,  after  samples  were  centrifuged  for  10 
min,  on  a  Hazen  platinum -cobalt  scale  with  a  Hellige  aqua  tester  model  611A.  Turbidity  was 
measured  with  a  Hach  Turbidimeter  model  2100 A  and  conductivity  was  measured  with  a 
conductivity  Bridge  meter  model  31. 

Data  Manipulation-  Average  monthly  TP,  TDP,  and  FPP  concentrations  were 
calculated.  These  averages  were  calculated  after  weighting  each  set  of  values  by  the 
corresponding  time  interval  between  sample  collection.  Phosphorus  loading  was  calculated  by 
summing  the  product  of  P  concentration  (mg/m3)  and  discharge  (m3/s)  for  all  sampling 
intervals.  Export  coefficients  were  calculated  by  dividing  the  total  mass  of  P  transported  by  the 
watershed  area.  Relationships  between  TP,  FPP,  and  TDP  concentrations  and  discharge  were 
analyzed  on  an  annual  and  seasonal  basis  with  regression  and  correlation  analysis.  To  reduce 
outliers,  these  analyses  were  done  after  the  data  were  transformed  to  natural  logarithms.  When 
the  data  were  grouped  seasonally  or  on  the  rising  vs.  falling  stage  of  the  discharge  hydrograph, 
analysis  of  covariance  (ANCOVA)  was  used  to  test  for  equality  of  the  regression  lines.  Linear 
regression  lines  were  compared  in  all  cases  because  (1)  not  all  relationships  were  described  by 
nonlinear  models,  (2)  nonlinear  and  linear  correlation  coefficients  were  very  similar.  Analysis 
of  covariance  was  also  used  to  compare  slopes  between  the  various  P  fractions  and  discharge. 
Degree  of  similarity  for  TP  loading  between  months  was  determined  by  cluster  analysis.  The 
regression,  ANCOVA,  and  cluster  analyses  were  performed  with  BMDP  statistical  software 
(Dixon  1981)  on  the  Amdahl  computer  model  580/5860  at  the  University  of  Alberta. 

Replicate  grab  samples  were  compared  with  paired  t- tests  (Snedecor  and  Cochran 
1980).  Data  collected  to  estimate  spatial  heterogeneity  as  well  as  to  test  for  diel  variation  were 
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analyzed  with  the  maximum  normed  residual  (MNR)  test  for  outliers.  Since  multiple 
comparisons  were  used  when  looking  at  diel  variation  and  spatial  heterogeneity,  significance 
was  accepted  at  P<0.01.  This  test  is  not  refined  enough  to  test  for  small  differences;  only  large 
differences  could  be  detected. 

Data  from  Two  Creek  were  used  for  most  illustrations  except  where  important 
differences  between  the  two  streams  warranted  added  illustrations.  Illustrations  for 
Sakwatamau  River  data  are  in  Appendix  C.  Raw  data  for  both  Two  Creek  and  Sakwatamau 
River  are  listed  in  Appendix  D. 


D.  Site  Description 

This  study  was  conducted  on  two  adjacent,  undisturbed  watersheds  (Two  Creek  and 
Sakwatamau  River)  located  approximately  230  km  northwest  of  Edmonton,  Alberta  (Fig.  2.1). 
At  the  sampling  sites.  Two  Creek  (54°18'N,  116020'W)  is  a  third-order  stream  draining  161 
km:  and  Sakwatamau  River  (54°24'N,  116°05'W)  is  a  fourth-order  stream  draining  a  281  km: 
watershed.  Two  Creek  and  Sakwatamau  River  drain  into  the  Athabasca  River  8  and  30  km 
downstream  from  the  sampling  sites,  respectively. 

Watershed  vegetation  is  mixed  boreal  forest;  the  dominant  tree  species  are  white  spruce 
( Picea  glauca),  lodgepole  pine  ( Pinus  contorta ),  aspen  poplar  ( Populus  tremuloides) ,  and 
balsam  poplar  ( Populus  balsamifera) .  Shrub  species  along  stream  banks  include  willows  ( Salix 
spp.),  wild  rose  ( Rosa  acicularis),  river  alder  ( Alnus  tenui folia),  and  various  berries.  Stream 
faunal  communities  are  represented  by  invertebrates  (primarily  Corixidae  and  Ephemeroptera; 
Appendix  A),  mountain  whitefish  ( Prosopium  williamsoni) ,  and  arctic  grayling  ( Thymallus 
arcticus ) . 

Bedrock  underlying  the  streams  is  sedimentary  (shale,  sandstone,  and  mudstone)  of  the 
Paskapoo  formation  (Paleocene)  overlain  by  3  to  5  m  of  glacial  till  (Knapik  and  Lindsay 
1983).  Watershed  soils  are  primarily  Gleysolic,  which  is  a  clay  or  sandy-clay  loam,  and  organic 
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soils  are  found  along  the  stream  beds  (Knapik  and  Lindsay  1983). 

Dominant  ions  in  the  streamwater  are  calcium  and  bicarbonate;  average  conductivity 
was  moderate  (185  and  174  wmhos/cm  for  Two  Creek  and  Sakwatamau  River,  respectively). 
Average  turbidity  of  streamwater  was  low  (2.5  and  2.6  NTU  for  Two  Creek  and  Sakwatamau 
River,  respectively)  although  average  colour  was  moderate  (69  and  47  mg/L  Pt  for  Two  Creek 
and  Sakwatamau  River,  respectively).  Total  nitrogen  to  total  P  ratios  were  moderate  (22-27). 
Average  pH  was  high  (8.1)  for  both  streams  and  water  temperature  ranged  between  0°  and 
20°C. 


E.  Results 

When  replicate  grab  samples  collected  at  the  same  time  and  location  were  compared, 
there  were  no  significant  differences  in  TP  concentrations  (16.6 ±1.2  mg/m3,  n=ll,  r  =  1 .82, 
^>0.05)  and  TDP  concentrations  (10.2±0.4  mg/m3,  «=11,  r=0.75,  ^>0.4)  in  both  streams. 
Similarly,  when  samples  collected  across  the  width  of  either  stream  at  approximately  the  same 
time  were  compared,  there  were  no  differences  in  FPP  or  TDP  concentrations  on  any  of  the 
dates  sampled  (P>0.05;  Fig.  2.2).  However,  there  were  significant  changes  in  discharge  for 
Two  Creek  on  May  11,  June  2,  and  August  30  (FCO.Ol).  Therefore,  in  the  two  study  streams, 
P  was  mixed  fairly  well  horizontally  such  that  even  when  discharge  was  so  low  it  was  not 
detectable  at  a  location  (e.g.  7-m  on  May  11  at  Two  Creek;  Fig.  2.2),  P  concentration  in  the 
water  column  was  uniform.  Thus,  the  grab  sample  technique  used  in  this  study  was  adequate  to 
estimate  P  concentration  in  streamwater. 

On  eight  of  the  nine  sampling  dates,  there  were  no  diurnal  changes  (P>0.01)  in 
discharge,  FPP,  and  TDP  concentrations  for  Two  Creek  and  Sakwatamau  River  (Fig.  2.3). 
However,  for  Sakwatamau  River  on  May  5-6,  there  was  a  significant  change  in  discharge 
(F<0.01)  due  to  increased  runoff  resulting  from  rain.  These  results  suggest  no  diel  patterns  in 
P  concentrations  in  Two  Creek  and  Sakwatamau  River. 
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During  the  sampling  year,  623  mm  of  precipitation  fell  on  the  watersheds,  most  (75%) 
during  June  and  July.  Relatively  little  precipitation  fell  during  spring  (9%)  and  autumn  (10%). 
This  pattern  of  rainfall  was  reflected  in  trends  of  stream  discharge  (Fig.  2.4).  Heavy  storms 
during  June  and  July  caused  dramatic  increases  in  discharge  accompanied  by  overland  flow. 
Consequently,  during  June  and  July,  82%  and  64%  of  annual  discharge  was  recorded  for  Two 
Creek  and  Sakwatamau  River,  respectively.  This  pattern  of  heavy  rainfall  during  early  summer 
is  not  unusual  for  the  region;  maximum  daily  discharge  at  the  mouth  of  Sakwatamau  River  was 
during  June  in  5  of  the  last  11  years.  Only  once  during  the  11 -yr  period  was  maximum  daily 
discharge  recorded  during  spring  runoff  (April  27,  1974).  Based  on  data  collected  by 
Environment  Canada  at  the  mouth  of  Sakawatamau  River,  total  annual  discharge  in  1983  (164 
000  dam3)  was  slightly  below  the  11 -yr  average  (179  145  dam3).  However,  the  maximum 
instantaneous  discharge  (the  largest  discharge  recorded  at  any  one  time  in  the  sampling  year)  in 
1983  was  exceeded  three  times  since  1973.  Since  overland  flow  accompanied  maximum 
instantaneous  discharge  at  the  study  site  in  1983,  overland  flow  probably  occurred  at  least  four 
times  during  the  last  11  years.  Thus,  materials  transported  during  overland  flow  could  have 
accumulated  for  several  years. 

Storm  Analysis  of  P  Concentrations-  Since  such  a  high  percentage  of  annual  rainfall  fell 
during  late  June  and  July,  discrete  storm  analysis  is  difficult  and  cumulative  effects  were  very 
important.  The  first  rainstorm  (May  6-8)  came  before  all  the  snow  in  the  watershed  had  melted 
and  consequently,  a  sustained  increase  was  recorded  in  discharge  and  TP  concentration  because 
of  increased  snowmelt.  The  next  storm  event  (June  17-21)  was  accompanied  by  overland  flow. 
During  this  storm,  TP  concentration  increased  by  851%  and  1491%  over  pre-storm 
concentrations  in  Two  Creek  and  Sakwatamau  River,  respectively;  87%  and  93%  of  the  increase 
was  in  the  FPP  fraction.  At  the  same  time,  CPM  transport  increased  by  1001%  in  Two  Creek 
and  by  529%  in  Sakwatamau  River  over  pre- storm  transport.  Peak  discharge  and  TDP 
concentration  were  measured  the  day  following  peak  rainfall.  However,  peak  TP  concentration 
was  noted  three  days  after  peak  discharge  and  consequently  TP  concentrations  were  higher  on 


•  I  q  J 


22 


the  falling  limb  of  the  storm  hydrograph  (hysteresis  effect  predicts  higher  P  concentrations  at 
similar  discharges  on  the  rising  limb  than  on  the  falling  limb  of  the  storm  hydrograph  (Gray 
1970)).  Storms  on  June  30  to  July  3  and  July  15  to  17  resulted  in  similar  trends  although 
increases  in  FPP  (443%  and  289%  in  the  two  storms,  respectively)  and  TDP  (134%,  126%) 
concentrations  were  not  so  dramatic. 

Rain  fell  on  most  days  in  July  but  heavy  rain  storms  on  July  25  to  27  and  again  on  July 
29  resulted  in  a  266%  increase  in  TP  concentration  with  93%  of  this  increase  in  the  FPP  fraction 
in  Two  Creek.  During  these  storms,  CPM  transport  increased  by  only  141%;  however,  the 
relationship  between  the  storm  hydrograph  and  P  concentration  showed  a  different  pattern 
from  previous  storms.  Total  P  concentrations  peaked  the  same  day  as  rainfall  while  stream 
discharge  and  TDP  concentrations  peaked  the  following  day,  although  changes  in  TDP 
concentrations  were  slight.  Consequently,  TP  concentrations  were  higher  on  the  rising  limb 
than  on  the  falling  limb  of  the  storm  hydrograph  (hysteresis  effect). 

Patterns  in  Phosphorus  Dynamics-  Average  TP,  TDP,  and  FPP  concentrations 
(time -weighted)  over  the  sampling  year  were  very  similar  in  Two  Creek  and  Sakwatamau  River 
(Table  2.1).  In  both  streams,  average  monthly  TP  and  FPP  concentrations  were  highest  in 
June  when  discharge  was  also  high.  In  Sakwatamau  River,  average  monthly  TDP  concentration 
was  highest  in  July  whereas  in  Two  Creek,  TDP  was  highest  in  April.  In  Sakwatamau  River, 
the  second  highest  average  monthly  TP,  TDP,  and  FPP  concentrations  were  all  in  April.  The 
reason  for  relatively  low  spring  FPP  concentration  in  Two  Creek  compared  to  Sakwatamau 
River  can  be  postulated.  Ice  on  the  Sakwatamau  River  streambed  melted  before  April  22  and  all 
ice  along  the  banks  was  gone  by  April  24.  However,  in  Two  Creek,  ice  melted  from  the  top  to 
the  bottom  and  ice  sitting  on  the  substrate  did  not  melt  until  May  3.  Fine  particulate  P 
transport  was  significantly  greater  in  Sakwatamau  River  than  in  Two  Creek  from  April  22  - 
May  3  (Wilcoxon  two-sample  rank  test,  «=14,  PC0.01).  This  suggests  that  there  was  no 
resuspension  of  FPP  from  the  streambed  in  Two  Creek  during  this  period  whereas  there  was 
resuspension  from  the  streambed  in  Sakwatamau  River.  In  both  streams,  TP,  TDP,  and  FPP 
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concentrations  were  lowest  in  October  when  discharge  was  also  low. 

Most  of  the  P  in  the  FPP  size  fraction  from  August  to  October  was  bound  to  materials 
in  the  5-10  um  size  range  for  both  Two  Creek  (50%)  and  Sakwatamau  River  (71%).  The  next 
most  important  fraction  in  FPP  was  the  0.45  to  5.0  um  fraction  in  Two  Creek  (33%)  and  in 
Sakwatamau  River  (12%),  especially  in  the  autumn.  There  was  relatively  little  P  in  the  10-100 
um  size  fractions  in  Two  Creek  (18%)  and  Sakwatamau  River  (17%)  from  August  to  October. 
From  April  to  June,  1984,  P  was  bound  to  all  fractions  of  FPP;  no  fraction  was  dominant 
(Figs.  6.3  and  6.4). 

Discharge  and  FPP  concentration  were  strongly  correlated  (r=0.81  and  r=0.74  for 
Two  Creek  (Fig.  2.5)  and  Sakwatamau  River(Fig  6.5),  respectively).  For  both  streams,  a 
second  order  polynomial  fits  the  FPP  -  discharge  (Q)  data  better  than  the  linear  model,  and 
when  the  rate  of  change  of  discharge  (R)  was  added  as  a  third  independent  variable,  the 
nonlinear  model  fits  the  FPP  data  even  better  (P<0.05).  In  both  streams,  when  the  FPP  on  Q 
and  R  relationship  was  compared  between  rising  and  falling  stage  of  the  hydrograph,  there  were 
no  significant  differences  (n=201,  7>=0.12  for  Two  Creek;  /j=165,  P=0.15  for  Sakwatamau 
River).  In  Two  Creek,  there  were  no  significant  seasonal  differences  in  the  relationship  of  FPP 
on  discharge  between  early  vs.  late  summer  storms  (n= 65,  F=0.08);  there  were  insufficient 
data  to  do  a  similar  analysis  for  Sakwatamau  River  data.  In  both  Two  Creek  (Fig.  2.6)  and 
Sakwatamau  River  (Fig.  6.6),  TDP  concentration  was  positively  correlated  to  discharge 
(r=0.68  and  r=0.66,  respectively)  although  the  slope  of  the  regression  of  TDP  on  Q  was 
significantly  lower  than  the  slope  of  the  regression  of  FPP  on  Q  (/,<0.01).  For  the  TDP  -  Q 
data,  a  linear  model  was  the  best  fit  (i^O.05).  Finally,  the  relationship  between  TP 
concentration  and  discharge  was  similar  to  the  relationship  between  FPP  and  discharge, 
although  more  of  the  variability  in  the  data  set  was  explained  with  the  TP  models  than  the  FPP 
models  in  both  Two  Creek  (Eq.l)  and  Sakwatamau  River  (Eq.  2). 


InTP  =  3.04  +  0.251nQ  +  0.09(lnQ)2  +  0.12R, 


(r=0.87,  n=  209) 


(1) 


•  .  . 


24 


InTP  =  2.91  +  0.261nQ  +  0.12(lnQ)2  -  0.09R,  (r=0.78,  n  =  166) 


(2) 


In  Two  Creek,  contrary  to  the  FPP  on  Q  and  R  results,  there  was  a  significant  difference 
between  regressions  of  TP  on  Q  and  R  based  on  early  (June  18  -  July  13;  Eq.  3)  vs.  late 
summer  storm  data  (July  14  *  August  4;  Eq.  4)  (P<0.001). 


InTP  =  2.36  +  0.541nQ  +  0.12(lnQ)2  -  0.33R,  (r=0.94,  n  =  33) 


(3) 


InTP  =  3.06  +  0.091nQ  +  0.14(lnQ)2  +  0.25R,  (r=0.90,  «  =  33) 


(4) 


The  mass  of  coarse  particulate  matter  transported  downstream  during  storms  was  large; 
48  533  and  5071  kg  dry  weight  was  transported  from  Two  Creek  and  Sakwatamau  River, 
respectively.  The  relationship  between  CPM  and  discharge  is  linear  reflecting  the  strong 
increase  in  total  CPM  transport  during  high  discharge  (r=0.88  and  r=0.79  for  Two  Creek  and 
Sakwatamau  River,  respectively;  Fig.  2.7).  The  seasonal  pattern  of  P  content  of  CPM  was 
consistent  between  streams;  CPM  transport  was  largest  in  June  (25  918  and  4365  kg  for  Two 
Creek  and  Sakwatamau  River,  respectively)  and  lowest  in  October  (3  kg  for  both  streams). 
Phosphorus  content  of  CPM  generally  increased  over  summer  (from  0.13%  to  0.21%  of  dry 
weight  in  Two  Creek  and  from  0.11  to  0.22%  of  dry  weight  in  Sakwatamau  River)  and  then 
declined  during  peak  needle  and  leaffall  in  autumn  (Table  2).  This  suggests  that  P  content  of 
CPM  increases  with  length  of  time  in  the  stream.  However,  99  and  98%  of  CPM  for  Two  Creek 
and  Sakwatamau  River,  respectively,  was  transported  during  storms  in  June  and  July  so  CPP 
loading  was  highest  in  summer  and  lowest  in  autumn  when  flows  were  low  (Tables  2.3  and 
2.4). 

During  the  sampling  period,  13.04  mg/m2  and  7.49  mg/m2  of  P  (including  TP  and 
CPP)  was  exported  from  Two  Creek  and  Sakwatamau  River,  respectively.  Both  streams  carried 
72%  of  the  annual  P  export  in  the  FPP  size  fraction.  Coarse  particulate  matter  was  a  relatively 
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insignificant  portion  of  the  P  load;  3%  and  0.3%  of  total  P  transport  for  Two  Creek  and 
Sakwatamau  River,  respectively.  In  Two  Creek,  P  transport  was  greatest  in  June  for  all  three 
size  fraction.  In  Sakwatamau  River,  this  general  pattern  also  held,  although  TDP  transport  was 
also  high  in  July.  For  all  size  fractions  in  both  streams,  P  loading  was  lowest  during  October. 
In  all  months  except  June  and  July,  TDP  transport  was  greater  than  FPP  transport  for  both 
streams.  To  determine  which  months  were  most  closely  related  in  terms  of  P  export,  a  cluster 
analysis  was  used;  TP  loading  in  the  months  of  June,  July,  and  August  were  most  closely 
related.  These  were  next  most  similar  to  May  and  then  September.  The  first  and  last  months  of 
the  study  (April  and  October)  were  the  least  related  to  any  other  month.  This  trend  held  true 
for  both  Two  Creek  and  Sakwatamau  River  (Fig.  6.8).  This  pattern  suggests  a  seasonal 
influence  in  TP  loading  reflecting  large  increases  in  discharge  during  the  summer  months. 

Phosphorus  export  was  strongly  pulsed  during  storm  events;  in  Two  Creek,  42%  of  TP 
was  exported  with  19%  of  the  water  in  2%  of  the  time  (Fig.  2.8A).  Export  of  FPP  was  more 
strongly  pulsed  than  was  export  of  TDP;  75%  of  FPP  was  transported  within  the  same  time 
period  as  47%  of  TDP  transport  (Fig.  2.8B).  To  determine  how  important  the  two  major 
storms  were  in  terms  of  pulsing,  the  mass  of  P  transported  during  storms  was  calculated.  In 
Two  Creek,  83%  of  the  TP  loading  in  June  (1097  Kg)  was  transported  during  the  first  storm 
(June  17-21)  and  this  storm  accounted  for  55%  of  TP  loading  for  the  entire  sampling  year.  In 
Two  Creek,  50%  of  the  TP  loading  in  July  (294  Kg)  was  transported  during  the  second  large 
storm  (July  24-30)  and  this  storm  accounted  for  14%  of  TP  transport  for  the  sampling  year. 
Consequently,  68%  of  TP  loading  was  transported  during  12  days  (6%  of  the  time);  76%  of 
FPP  loading  was  transported  during  the  two  storms.  Thus,  P  export,  primarily  in  the  FPP 
fraction,  was  strongly  pulsed  during  storm  events. 
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F.  Discussion 

Although  de  March  (1975)  reported  diel  changes  in  P  concentration  in  arctic  streams,  I 
found  no  diel  patterns  in  P  concentration  in  Two  Creek  and  Sakwatamau  River.  Meyer  and 
Likens  (1979)  also  observed  no  diel  trends  in  TDP  concentrations.  O'Connell  and  Davies 
(1977)  reported  differences  in  fine  particulate  detritus  density  over  a  24 -h  period  for  mountain 
headwater  streams  in  Alberta,  possibly  due  to  changes  in  the  rate  of  detrital  processing  by 
stream  invertebrates.  Kaplan  and  Bott  (1982)  postulated  several  reasons  for  observed  diel 
fluctuation  in  dissolved  organic  carbon:  diel  changes  in  inputs  of  allochthonous  litter, 
autochthonous  production  and  processing,  or  fluctuation  in  inputs  of  groundwater.  If  these 
processes  can  also  explain  diel  P  patterns,  I  postulate  that  other  physical  processes  such  as  flow 
were  dominant  over  detrital  processing  by  invertebrates  in  Two  Creek  and  Sakwatamau  River. 
Diel  patterns  may  be  more  prevalent  in  smaller  streams  where  discharge  is  low  and  thus  not 
likely  to  be  a  dominant  influence  or  in  streams  that  are  influenced  by  other  diel  patterns  such  as 
beaver  activity  or  point -source  effluents. 

Although  Meyer  and  Likens  (1979)  reported  no  relationship  between  TDP  and 
discharge  in  the  forested  Hubbard  Brook  streams,  there  was  a  positive  relationship  between 
TDP  and  discharge  in  Two  Creek  and  Sakwatamau  River.  However,  my  results  concur  with 
results  obtained  in  grazed  watersheds  in  Dartmoor,  England  (Rigler  1979).  In  the  study 
streams,  positive  correlations  between  TDP  concentration  and  discharge  could  be  a  result  of 
several  factors:  1)  higher  TDP  concentrations  at  higher  discharges  during  spring  caused  by 
resuspension  of  sediment,  that  had  accumulated  during  low  flows  the  previous  autumn,  from  a 
thawing  streambed,  2)  high  TDP  concentration  due  to  increased  allochthonous  inputs  with 
overland  flow  (at  high  discharge),  or  3)  decreasing  TDP  concentration  during  declining  flows 
in  autumn,  perhaps  due  to  increased  processing  and  thus  conversion  into  FPM  during  low 
flows.  These  explanations  would  suggest  a  seasonal  relationship  between  TDP  and  discharge. 
This  last  hypothetical  process  concerning  increased  processing  at  low  flows  is  supported  by 
research  on  P  spiralling  which  suggests  that  processing  within  stream  biotic  compartments  is 
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most  important  at  low  flows  (Newbold  et  al.  1983).  Meyer  and  Likens  (1979)  attribute  the 
absence  of  a  relationship  between  TDP  and  discharge  to  very  low  concentrations  of  TDP 
resulting  from  biotic  (e.g.  periphyton)  and  abiotic  (e.g.  sediments)  P  sinks.  However,  the 
Dartmoor  streams  have  an  organic -rich  soil  horizon  which  could  act  as  a  P  sink  which  would 
suggest  no  relationship  between  TDP  and  discharge,  contrary  to  the  observed  positive 
relationship  between  TDP  and  discharge.  In  Two  Creek  and  Sakwatamau  River,  periphyton 
only  grew  for  a  2-wk  period  in  early  fall  and  there  were  no  bryophytes  at  any  time  so  a  biotic 
sink  is  an  unlikely  explanation  in  these  streams.  During  low  flow,  stream  bacteria  may  have 
been  an  important  P  sink,  however,  bacteria  were  not  quantified  in  this  study.  A  possible 
explanation  for  a  positive  relationship  between  TDP  and  discharge  in  the  study  streams  is 
increase  in  groundwater  inputs  during  storm  events;  groundwater  is  usually  high  in  TDP.  A 
positive  relationship  may  exist  between  TDP  and  discharge  for  all  streams  overlying 
sedimentary  rock. 

Although  both  TDP  and  FPP  concentrations  increased  with  increasing  discharge  in  the 
two  study  streams,  increases  in  FPP  were  stronger  than  increases  in  TDP  (/><0.01).  These 
results  concur  with  earlier  studies  which  reported  a  stronger  response  of  FPP  concentration  to 
discharge  during  summer  than  winter  (Meyer  and  Likens  1979,  Rigler  1979).  As  well,  these  two 
earlier  studies  did  not  report  seasonal  differences  except  between  summer  and  winter  or  between 
rising  and  falling  stage.  In  contrast,  there  were  no  significant  differences  between  P 
concentrations  on  the  rising  vs.  falling  stage  of  the  discharge  hydrograph  in  Two  Creek  and 
Sakwatamau  River.  In  addition,  in  the  central  Alberta  streams,  storms  during  the  early  part  of 
the  summer  did  not  elicit  the  hysteresis  effect.  These  two  unexpected  results  suggest  that  a  large 
P  pool  was  still  contributing  to  P  in  the  flowing  water  while  discharge  was  decreasing.  A  large  P 
pool  could  have  accumulated  during  autumn  and  winter  and  remained  in  the  stream  during  the 
relatively  small  spring  runoff  compared  to  other  areas.  As  well,  P  may  have  accumulated  since 
spring  runoff.  In  the  later  part  of  the  summer,  hysteresis  effects  were  apparent  in  Two  Creek 
and  Sakwatamau  River  during  storms.  This  late  appearance  of  the  hysteresis  effect  suggests 
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that  during  the  early  summer  storms,  FPM  that  had  accumulated  over  a  considerable  time 
period  was  washed  out  of  the  streambeds  whereas  increases  in  later  summer  were  due  to  P 
which  had  accumulated  over  a  much  shorter  time  period,  i.e.,  since  the  last  storm,  as  was 
suggested  by  Rigler  (1979).  Thus,  all  these  detailed  stream  studies  indicate  that  FPP  transport 
is  dependent  on  discharge  as  well  as  length  of  time  since  the  last  storm.  The  dependency  of  FPP 
transport  on  discharge  is  further  supported  by  FPP  fractionation;  P  was  transported  in  all 
fractions  at  high  flows  when  coarser  sediments  were  suspended  and  in  smaller  fractions  during 
low  flows  when  only  finer  sediments  were  transported . 

In  Two  Creek  and  Sakwatamau  River,  P  content  of  CPM  increased  during  summer  and 
peaked  during  autumn.  This  pattern  was  also  observed  on  hardwood  watersheds  at  Hubbard 
Brook  (Meyer  and  Likens  1979).  However,  there  was  little  seasonal  variation  in  nutrient 
content  of  coniferous  and  deciduous  fractions  from  two  streams  in  Ontario  (Dance  et  al.  1979). 
Since  the  watersheds  in  this  study  are  dominated  by  coniferous  trees,  leaffall  was  slight  but 
large  amounts  of  spruce  needles  were  transported  in  stream  water.  In  the  study  streams,  the 
decline  in  P  content  of  CPM  at  peak  leaffall  could  be  caused  by  retention  of  nutrients  by  trees 
when  they  shed  their  needles  or  leaves.  Even  though  Dance  et  al.  (1979)  reported  lower  P 
content  in  coniferous  matter  than  deciduous  matter,  the  P  content  of  CPM  from  Two  Creek 
and  Sakwatamau  River  (0.10%  -  0.22%)  was  higher  than  for  a  deciduous  forest  stream 
(0.079%)  at  Hubbard  Brook  (Gosz  et  al.  1972).  Geographic  variation  could  explain  this 
discrepancy.  Even  though  the  P  content  of  CPM  was  highest  in  autumn  in  the  two  study 
streams,  P  transport  of  CPM  was  greatest  in  summer  because  of  the  large  mass  of  CPM 
transported  during  storms.  While  CPP  was  transported  primarily  in  autumn  at  Hubbard  Brook, 
this  autumn  peak  can  be  explained  by  leaffall  from  the  deciduous  watershed  and  increased 
flows  in  autumn.  The  two  different  trends  suggest  that  CPP  export  is  dependent  upon  type  of 
watershed  vegetation  and  regional  pattern  of  discharge. 

Export  coefficients  from  Two  Creek  and  Sakwatamau  River  were  13.0  mg/m2  and  7.5 
mg/m2,  respectively.  Export  coefficients  in  this  study  are  similar  to  coefficients  developed  for 
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forested  watersheds  overlying  sedimentary  bedrock  in  southern  Ontario  (11.7  mg/m2;  Dillon 
and  Kirchner  1975).  While  Dillon  and  Kirchner's  (1975)  coefficient  is  relatively  high,  partly 
because  of  large  inputs  in  autumn  from  deciduous  trees,  their  estimate  is  probably  an 
underestimate  because  of  inadequate  sampling  frequency.  As  expected,  export  from  the  two 
undisturbed  study  streams  was  less  than  export  from  grazed  watersheds  in  Dartmoor  (19.9 
mg/m2;  Rigler  1979)  and  Ireland  (25.1  mg/m2;  Stevens  and  Smith  1978).  Since  Dillon  and 
Kirchner's  (1975)  study  reported  less  P  export  from  igneous  than  sedimentary  watersheds,  it  is 
not  surprising  that  P  export  from  the  Hubbard  Brook  stream  (7.3  mg/m2;  Meyer  and  Likens 
1979)  was  less  than  from  Two  Creek  and  Sakwatamau  River.  One  major  difference  between 
this  study  and  the  Hubbard  Brook  study  was  that  CPP  accounted  for  19%  of  P  export  from  the 
watershed  in  Hubbard  Brook  whereas  CPP  accounted  for  3%  and  0.3%  of  the  export  from  Two 
Creek  and  Sakwatamau  River,  respectively.  One  major  factor  that  could  explain  this  difference 
is  autumn  leaf  fall  from  deciduous  trees  in  Hubbard  Brook.  In  contrast.  Two  Creek  and 
Sakwatamau  River  watersheds  are  forested  primarily  by  spruce  trees  and  there  is  no  autumn 
leaf  fall.  This  study  supports  the  general  hypothesis  that  P  export  is  predictable  from  bedrock 
geology  and  land  use  in  the  watershed. 

However,  there  is  no  apparent  reason  why  the  Two  Creek  export  coefficient  is  69% 
higher  than  the  Sakwatamau  River  export  coefficient.  Average  monthly  P  concentrations  were 
similar  for  the  two  streams  and  P  mass  exported  by  the  two  streams  were  also  very  similar; 
Sakwatamau  River  exported  3%  more  P  than  Two  Creek  based  on  absolute  mass.  The  relatively 
large  watershed  size  results  in  a  low  export  coefficient  for  Sakwatamau  River  as  compared  to 
Two  Creek.  This  discrepancy  suggests  that  the  actual  size  of  the  Sakwatamau  River  watershed 
may  be  smaller  than  that  calculated  from  surficial  topography.  As  well,  groundwater  may  alter 
P  concentrations  or  the  direction  of  subsurface  flow.  These  data  suggest  that  export  coefficients 
based  on  watershed  size  may  not  be  appropriate  for  streams  on  glacial  till  over  sedimentary- 
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In  the  study  streams,  P  transport  was  strongly  pulsed  during  the  first  part  of  the 
summer  and  storm  events  were  the  dominant  factor  influencing  P  transport.  Other  studies 
report  strong  P  pulsing,  although  sources  and  timing  of  the  pulses  vary.  For  example,  in  Bear 
Brook,  a  forested  stream  in  New  Hampshire,  67 %  of  the  annual  exports  were  transported 
during  10  days,  primarily  during  spring  and  autumn  (Meyer  and  Likens  1979).  In  two 
contour -farmed  watersheds,  50%  of  P  exports  were  transported  during  snowmelt  runoff  in 
March  (Alberts  et  al.  1978).  In  a  rural  stream  in  New  York,  75%  of  the  P  was  exported  during 
highest  flows  in  10%  of  the  time  (Johnson  et  al.  1976).  A  large  portion  of  stream  P  inputs  can 
be  exported  within  the  first  few  days  of  spring  runoff  for  Char  Lake  (de  March  1975)  and  few 
hours  for  the  Saqvaqjuak  region  (H.E.  Welch,  Freshwater  Institute,  pers.comm.)  in  the 
Canadian  Arctic.  Low  winter  snowpack  in  central  Alberta  compared  to  southern  Ontario, 
Quebec,  and  the  New  England  states  as  well  as  a  long  period  of  snowmelt  (e.g.  a  Chinook 
followed  by  a  snowstorm  in  April  1983)  could  account  for  the  negligible  importance  of  spring 
runoff  compared  to  other  studies.  Regional  differences  in  patterns  of  storm  events  could 
account  for  some  differences  in  timing  of  P  transport.  For  example,  in  New  Hampshire,  major 
storms  are  more  likely  to  occur  in  spring  and  autumn  (Meyer  and  Likens  1979)  whereas  in 
central  Alberta,  major  storms  are  more  likely  to  occur  in  early  summer.  Thus,  the  timing  of  P 
pulsing  may  be  predictable  on  a  regional  scale  based  on  climatic  regimes  and  this  information 
can  be  used  to  plan  monitoring  programs  for  sampling  fluvial  P.  For  example,  this  study 
indicates  that  a  major  portion  of  the  P  loading  takes  place  during  storm  events  and  that 
sampling  frequency  should  be  on  the  order  of  every  few  hours  during  heavy  storms  in  central 
Alberta.  As  well,  in  this  region,  daily  sampling  is  adequate  during  the  low  spring  runoff  and 
bi-monthly  sampling  during  low  flow  in  the  fall.  However,  the  efficiency  of  this  suggested 
sampling  regime  is  strongly  dependent  on  climatic  patterns  and  year-to-year  variation  may  be 
significant  (Meyer  and  Likens  1979).  When  TP  loading  coefficients  were  recalculated  based  on 
a  weekly  sampling  frequency,  the  coefficients  for  Two  Creek  ranged  between  an  underestimate 
of  54%  to  an  overestimate  of  10%  (mean  underestimate  of  15%),  depending  on  which  day  of  the 
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week  sampling  was  done.  Coefficients  based  on  weekly  sampling  for  Sakwatamau  River  were  all 
underestimates  ranging  between  10  and  49%  (mean  34%).  This  underestimation  indicates  the 
inadequacy  of  weekly  sampling  regimes.  As  well,  the  discrepancy  between  the  coefficients  for 
the  two  streams  based  on  weekly  samples  indicates  the  unpredicatability  of  weekly  sampling 
regimes. 

The  degree  and  causes  of  pulsing  have  important  implications  for  the  biotic  community 
within  the  stream.  If  pulsing  is  caused  by  large  increases  in  discharge,  as  in  this  study,  P  cycling 
distance  would  be  long  at  high  discharges  since  distance  travelled  in  one  cycle  is  the  product  of 
cycling  time  and  downstream  velocity  (Newbold  et  al.  1983).  Although  P  uptake  rate  does 
increase  at  higher  flows  because  of  turbulent  contact  of  the  water  column  with  the  substrate, 
this  increase  would  likely  be  overshadowed  in  the  two  study  streams  by  the  rapid  downstream 
transport  of  FPM  and  CPM  at  high  flows  and  the  large  increase  in  water  depth.  Estimations  of 
biologically  available  P  for  streams  in  Quebec  indicated  low  availability  of  the  suspended 
particulate  fraction  (Peters  1981)  which  is  the  largest  fraction  transported  during  high  storm 
flow  in  Two  Creek  and  Sakwatamau  River  (see  Chapter  5).  During  low  flows  in  all  streams,  P 
cycling  length  would  be  shorter  and  a  P  atom  would  be  retained  longer  in  the  stream  within 
both  biotic  and  abiotic  compartments.  Thus,  while  estimating  mass  of  P  transported 
downstream  is  undoubtedly  important  to  downstream  water  bodies,  mass  transported  does  not 
indicate  availability  of  P  to  the  stream  community  itself.  Rather,  P  uptake  and  cycling 
(assuming  some  degree  of  relatedness  to  availability)  is  dependent  upon  a  number  of  factors 
including  flow  rates,  standing  stock  of  particulate  matter,  and  standing  stock  of  exchangeable 
P. 

This  study  concurs  with  earlier  predictions  that  abiotic  factors,  such  as  watershed  land 
use,  bedrock  geology,  and  climatic  regimes,  all  influence  P  export  patterns.  However,  this  study 
suggests  that  P  export  coefficients  based  on  watershed  area  may  not  be  an  appropriate  export 
expression  for  forested  streams  overlying  sedimentary  bedrock.  Alberta  streams  would  be  ideal 
for  many  studies  including:  (1)  a  P  budget  for  headwater  streams  on  coniferous  watersheds 
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where  the  CPP  inputs  and  outputs  are  low  relative  to  deciduous  watersheds  and  where 
groundwater  may  be  a  significant  input  to  the  system,  (2)  availability  of  P  during  high  vs.  low 
flows,  (3)  possible  retention  or  accumulation  of  P  in  stream  sediments,  and  (4)  the  influence 
of  watershed  size  on  P  export  in  streams  overlying  glacial  till. 
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Table  2.1.  Time-weighted  average  total  phosphorus  (TP),  fine  particulate  phosphorus 
(FPP),  and  total  dissolved  phosphorus  (TDP)  concentrations  (mg/m3)  for  Two  Creek 
and  Sakwatamau  River. 


Month 

TWO  CREEK 

SAKWATAMAU  RIVER 

TP 

TDP 

FPP 

TP 

TDP 

FPP 

April 

28.4 

17.5 

10.9 

37.0 

21.7 

15.2 

May 

22.0 

11.6 

10.3 

22.7 

12.7 

10.1 

June 

43.3 

11.6 

32.2 

48.3 

11.0 

38.2 

July 

35.0 

14.7 

20.3 

32.2 

23.3 

13.1 

August 

20.5 

11.2 

9.4 

18.0 

10.7 

7.4 

September 

15.8 

9.3 

6.6 

14.6 

9.7 

5.0 

October 

13.5 

8.6 

4.9 

13.6 

9.2 

4.4 

Average 

25.5 

11.7 

14.0 

25.8 

13.1 

13.1 
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Table  2.2.  Percent  phosphorus  (P)  content  of  coarse  particulate  matter  (CPM)  for  Two 
Creek  and  Sakwatamau  River  for  5  dates  in  1983. 


DATE 

TWO  CREEK 

SAKWATAMAU  RIVER 

June  17 

0.13 

0.11 

August  1 

0.14 

0.16 

September  6 

0.16 

0.16 

September  22 

0.21 

0.22 

October  10 

0.10 

0.12 
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Table  2.3.  Mass  of  total  phosphorus  (TP),  total  dissolved  phosphorus  (TDP),  fine 
particulate  phosphorus  (FPP),  and  coarse  particulate  phosphorus  ( CP P)  in  kg  transported 
past  Two  Creek  study  site  during  each  month  of  the  study  and  over  the  whole  study  period 
expressed  in  total  mass  and  mass  per  unit  area  of  watershed  (Export  Coefficient  mg/m2). 
Note  that  TP  does  not  include  CPP . 


Month  TP  TDP  FPP  CPP 


April 

11.1 

May 

59.5 

June 

1326.0 

July 

584.4 

August 

37.5 

September 

11.3 

October 

5.4 

Total  Transport 

2035.2 

Export  Coefficient 

12.64 

6.7 

4.4 

29.8 

29.7 

249.5 

1076.5 

33.69 

201.2 

383.2 

30.43 

19.5 

18.0 

0.09 

6.6 

4.7 

0.01 

3.4 

2.0 

0.00 

516.7 

1518.5 

64.23 

3.21 

9.43 

0.40 

i  t'  it  5JC  1 
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Table  2.4.  Mass  of  total  phosphorus  (TP),  total  dissolved  phosphorus  (TDP),  fine 
particulate  phosphorus  (FPP),  and  coarse  particulate  phosphorus  (CPP)  in  kg  transported 
past  Sakwatamau  River  study  site  during  each  month  of  the  study  and  over  the  whole  study 
period  expressed  in  total  mass  and  mass  per  unit  area  of  watershed  (Export  Coefficient 
mg/m2).  Note  that  TP  does  not  include  CPP. 


Month 

TP 

TDP 

FPP 

CPP 

April 

53.1 

29.9 

23.2 

May 

86.7 

44.7 

42.0 

June 

1265.8 

181.2 

1084.6 

4.80 

July 

556.8 

191.2 

365.6 

0.83 

August 

76.4 

42.4 

33.9 

0.13 

September 

37.7 

25.2 

12.6 

0.02 

October 

22.8 

15.4 

7.4 

0.00 

Total  Transport 

2099.3 

530.0 

1569.3 

5.78 

Export  Coefficient 

7.47 

1.87 

5.58 

0.02 
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Fig.  2.1.  Map  of  study  area. 
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Fig.  2.2.  Discharge,  fine  particulate  phosphorus  (FPP).  and  total  dissolved  phosphorus 
(TDP)  concentrations  at  2-m  intervals  across  the  width  of  Two  Creek. 
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Fig.  2.3.  Diel  patterns  in  discharge  (Q).  total  phosphorus  (TP),  and  total  dissolved 
phosphorus  (TDP)  concentrations  for  Two  Creek  (T)  and  Sakwatamau  River  (S).  Note 
that  discharge  scale  changes.. 
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Fig.  2.5.  Relationship  between  fine  particulate  phorphorus  (FPP)  concentration  and 
discharge  (Q)  in  Two  Creek  (n  =  218).  R  is  rate  of  change  of  discharge. 
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Fig.  2.6.  Relationship  between  total  dissolved  phosphorus  (TDP)  concentration  and 
discharge  (Q)  in  Two  Creek. 
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Fig.  2.7.  Relationship  between  coarse  particulate  matter  (CPM)  transport  and  discharge 
(Q)  in  Two  Creek. 
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Fig.  2.8.  Flow  duration  curves  for  Two  Creek  indicating:  A.  the  %  total  time,  %  water 
output,  and  %  total  phosphorus  (TP)  transport  measured  at  or  above  the  specified 
discharge,  and  B.  the  %  TP,  %  fine  particulate  phosphorus  (FPP),  and  %  total  dissolved 
phosphorus  (TDP)  transport  measured  at  or  above  the  specified  discharge  in  Two  Creek. 
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III.  General  Discussion 


In  this  study,  epistemological  difficulties  involved  in  relating  patterns  with  processes 
(Hart  1983)  are  apparent.  For  example,  I  have  demonstrated  that  there  are  seasonal  patterns  in 
phosphorus  (P)  transport,  a  potential  limiting  energy  source  for  stream  biotic  communities. 
While  I  can  speculate  on  how  this  influences  functional  feeding  group  dynamics,  the 
predictions  would  involve  too  many  assumptions  to  be  useful.  These  assumptions  include:  (1) 
P  is  in  limiting  supply  to  at  least  the  dominant  species  or  functional  groups,  (2)  that  P  is 
limiting  at  all  times,  and  (3)  it  is  possible  to  isolate  the  influence  of  P  from  other  physical 
factors  such  as  seasonal  temperature  change.  An  enormous  effort  would  be  required  to  support 
these  assumptions.  Thus,  the  link  between  stream  community  processes  and  P  dynamics  is 
largely  speculative.  Stream  community  ecology  remains  a  young  science  although  the  prevailing 
view  is  that  abiotic  variables  (such  as  discharge)  override  biotic  variables  (such  as  invertebrate 
detritus  processing)  in  determining  stream  community  organization  (Barnes  and  Minshall 
1983). 

The  use  of  radioactive  tracers  has  been  important  in  discerning  the  role  of  P  within 
biotic  and  abiotic  compartments  within  streams  (Newbold  et  al.  1983).  The  cycling  (spiralling) 
of  P  includes  three  processes:  (1)  downstream  transport,  (2)  turnover  length  within  the  water 
column,  and  (3)  consumer  turnover  length.  However,  work  needs  to  be  done  on  metabolic  P 
requirements  of  stream  organisms  in  relation  to  P  availability  and  natural  variation  in 
metabolic  requirements.  Also,  the  interactions  of  P  in  flowing  water  with  the  sediments  and 
groundwater  has  largely  been  ignored  in  stream  P  studies. 

However,  the  detection  of  patterns  is  a  necessary  step  in  the  modelling  of  stream 
ecosystems.  This  study  represents  the  first  detailed  study  of  P  dynamics  on  coniferous 
watersheds;  all  previous  work  has  been  done  on  hardwood  watersheds  (e.g.  Meyer  and  Likens 
1979,  Newbold  et  al.  1983)  or  on  disturbed  watersheds  (e.g.  Rigler  1979).  Major  differences  in 
this  study  compared  with  previous  work  include  the  larger  proportion  of  P  being  transported  as 
dissolved  P  rather  than  fine  particulate  P  at  low  flows  (e.g.  63%  of  total  P  transport  in  October 
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was  in  the  dissolved  fraction)  and  the  relatively  limited  importance  of  coarse  particulate  P 
transport  (3%  of  total  P  transport).  Whether  these  differences  can  be  attributed  to  differences 
in  watershed  vegetation  is  unknown.  As  well,  regional  climatic  variation  can  alter  timing  of 
storms,  leaffall,  and  other  factors  that  have  a  strong  influence  on  P  export  dynamics  (Johnson 
et  al.  1976,  Meyer  and  Likens  1979).  The  models  (Eq.l  and  Eq.  2  whereby  Q  is  discharge  and  R 
is  rate  of  change  of  discharge)  I  constructed  for  two  streams  draining  coniferous  watersheds 

Two  Creek: 

InTP  =  3.04  +  0.251nQ  +  0.09(lnQ)2  +  0.12R  (n  =  209,  r=0.87)  (1) 

Sakwatamau  River : 

InTP  =  2.91  +  0.261nQ  +  0.12(lnQ)2  -  0.09R  (n  =  166,  r=0.78)  (2) 

can  be  refined  by  testing  data  obtained  from  other  coniferous  watersheds  that  differ  in  climatic 
patterns . 

The  detection  of  patterns  is  also  important  on  a  regional  scale  for  use  in  management 
programs.  For  example,  this  study  indicates  that  major  P  loading  occurs  during  storm  events 
which  suggests  that  sampling  frequency  should  be  on  the  order  of  every  few  hours  during  heavy 
summer  storms  in  central  Alberta.  As  well,  daily  sampling  would  be  adequate  during  spring 
runoff  and  bi-monthly  sampling  would  be  adequate  during  low  flows  in  the  fall.  However,  the 
efficiency  of  this  suggested  sampling  regime  is  strongly  dependent  on  climatic  patterns  and 
year-to-year  variation  may  be  significant  (Meyer  and  Likens  1979). 

Current  stream  ecosystem  models  are  descriptions  of  local  patterns  and  ignore  ecosytem 
stability  (Webster  et  al.  1983).  Stability  can  be  defined  in  terms  of  resistance  (ability  of 
ecosystem  to  exhibit  minimal  response  to  disturbance)  and  resilience  (rapidity  with  which  an 
ecosystem  responds  once  it  has  been  displaced  by  a  disturbance)  (Webster  et  al.  1983). 
Eventually,  stream  ecosystem  models  should  be  developed  to  predict  biotic  response  to 
stochastic  disturbances  within  the  ecosystem  based  on  the  magnitude  and  duration  of  the 
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disturbance.  The  rapid  return  to  pre- storm  P  export  patterns  exhibited  by  Two  Creek  and 
Sakwatamau  River  following  major  storms  would  suggest  a  resilient  ecosystem.  Even  though 
there  was  overland  flow  and  physical  alteration  of  the  stream  channels  (uprooted  trees  and 
streambank  erosion;  i.e.  strong  disturbance),  these  events  only  represent  a  short-term 
disturbance  (2  weeks).  Two  Creek  and  Sakwatamau  River  may  not  be  so  resilient  to  long-term 
disturbances.  Resistance  was  low  in  the  watersheds  since  stream  discharge  responded  strongly  to 
rain  events.  However,  evaluation  of  resistance  and  resilience  is  subjective  at  this  point  for 
stream  ecosystems. 

I  conclude  that  this  study  represents  an  important  advance  in  the  modeling  of  stream  P 
dynamics.  It  is  the  first  model  based  on  detailed  data  from  watersheds  with  the  following 
characteristics:  (1)  coniferous  trees,  (2)  3-5  m  of  glacial  till  overlying  sedimentary  bedrock, 
and  (3)  climatic  regime  with  strong  summer  storms.  The  direct  effect  of  each  characteristic  is 
unknown  although  the  combined  effect  is  important.  Phosphorus  transport  was  pulsed  during 
summer  storms;  68%  of  total  P  export  was  transported  during  12  days  for  Two  Creek.  While 
other  studies  on  both  forested  and  agricultural  watersheds  report  stronger  pulsing  (e.g.  Johnson 
et  al.  1976),  most  studies  report  pulsing  primarily  during  spring  and  autumn  (e.g.  Meyer  and 
Likens  1979).  The  P  pulsing  during  early  summer  storms  suggests  a  depletion  of  P  pools  during 
storm  events  and  an  accumulation  of  P  on  the  sediments  in  autumn.  The  potential  for  retention 
of  P  during  autumn  and  winter  in  Two  Creek  and  Sakwatamau  River  sediments  suggests  an 
energy  source  for  organisms  over -wintering  under  the  ice  and  in  the  sediments.  Whether  P 
retention  occurs  and  whether  P  retained  is  available  to  biotic  compartments  remain  to  be  tested. 
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IV.  APPENDIX  A:  Taxonomic  list  of  invertebrates  collected  from  Two  Creek  and  Sakwatamau 


River  using  dipnets. 


ORDER 

FAMILY 

GENERA 

Ephemeroptera 

Heptageniidae 

Heptagenia 

Cinygmula 

Ephemerellidae 

Ephemerella 

Siphionuridae 

Leptophlebiidae 

Ameletus 

Plecoptera 

Taeniopterygidae 

Doddsia 

Nemouridae 

Zapada 

Chloroperlidae 

U  taper  la 

Perlodidae 

I  so  per  la 

Perlidae 

Claasenia 

Acroneuria 

Pteronarcyidae 

Capniidae 

Pteronarcys 

Trichoptera 

Hydropsychidae 

Arcopsyche 

Brachycentridae 

Brachycenlrus 

Limnephilidae 

Limnephilus 

Hydalopylax 

Rhyacophilidae 

Rhyacophila 

Lepidostomatidae 

Lepidostoma 

T helio  psyche 

Diptera 

Tabanidae 

Chironomidae 

Tab anus 

Hemiptera 

Coleoptera 

Corixidae 
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V.  APPENDIX  B:  Longitudinal  transport  of  phosphorus 

The  study  of  material  transport  past  one  stream  site  yields  valuable  information  about 
temporal  changes  in  transport  and  the  influence  of  physical  characteristics  in  the  watershed  on 
transport  dynamics.  However,  the  measurement  of  changes  in  concentrations  and  relative 
importance  of  various  factors  as  materials  are  transported  downstream,  can  yield  information 
on  processes  occurring  within  the  stream.  Vannote  et  al.  (1980)  elaborated  on  these  ideas  and 
developed  a  model,  called  the  river  continuum  concept,  to  describe  the  structure  and  function 
of  communities  along  a  river  system.  Minshall  (1978)  reported  a  switch  in  dependence  from 
allochthonous  to  autochthonous  sources  of  energy  in  higher -order  woodland  streams  by 
measuring  photosynthesis  and  respiration  rates  as  well  as  inputs  and  exports  of  organic  matter 
which  supports  the  river  continuum  model.  A  major  implication  is  that  downstream 
communities  capitalize  on  upstream  processing  inefficiencies  and  that  both  upstream 
inefficiencies  and  downstream  adjustments  are  predictable.  General  predictions  of  the  river 
continuum  model  were  congruent  with  results  obtained  from  analyses  of  longitudinal  and 
seasonal  shifts  in  functional  group  abundance  with  concurrent  changes  in  food  availability  in 
Oregon  streams  (Hawkins  and  Sedell  1981)  and  in  an  Albertan  river  system  (Culp  and  Davies 
1982). 

The  measurement  of  food  availability  includes  the  measurement  of  materials  being 
transported  within  different  size  fractions.  Food  availability  is  then  usually  related  to  the 
density  of  animals  at  each  site  most  likely  to  use  that  size  fraction  as  a  food  source  (i.e. 
functional  group);  for  example,  shredders  vs.  coarse  particulate  organic  matter.  Thus,  for  each 
site  along  a  river,  the  relationship  between  potential  food  resource  and  invertebrate  density  is 
known  and  can  be  tested  against  predictions  from  the  continuum  hypothesis.  The  biggest 
problem  with  this  methodology  is  the  lack  of  knowledge  about  cause  and  effect.  For  example,  is 
a  shift  in  abundance  of  a  functional  group  a  response  to  a  shift  in  the  size  fraction  of  materials 
being  transported,  or  is  the  distribution  of  animals  a  result  of  changing  substrate  types?  Is  the 
shift  in  size  of  materials  being  transported  a  reflection  of  different  processing  rates  caused  by  a 
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change  in  the  abundance  of  a  functional  group?  Macroinvertebrates  are  an  important 
component  in  the  longitudinal  linkage  because  their  influence  on  the  quantity  and  timing  of 
fine  particulate  matter  transport  (Webster  1983).  The  obvious  approach  to  overcome  this 
problem  is  well -controlled  field  experiments  that  alter  either  food  availability  or  abundance  of  a 
functional  group  and  measure  the  response  within  other  ecosystem  components.  The  focus 
would  be  more  on  process  than  pattern. 

On  August  4  and  October  4  1983,  I  measured  phosphorus  (P)  partitioning  at  9  sites 
(plus  3  seeps  which  are  springs  or  small  streams  draining  into  the  main  channel)  along  a  3390  m 
section  of  Two  Creek  to  see  if  it  was  possible  to  detect  changes  in  the  partitioning  of  P  en  route 
downstream  (Fig.  5.1).  Phosphorus  was  separated  into  two  size  fractions:  total  dissolved 
phosphorus  (TDP,  <0.45um)  and  fine  particulate  phosphorus  (FPP,  0.45um  -  1mm).  The  P 
energy  pathway  was  chosen  to  test  for  changes  in  food  availability  for  two  reasons:  (1) 
measurement  of  benthic  biomass  and  transport  of  material  is  not  necessarily  an  adequate 
estimate  of  food  availability  because  food  quality  varies,  and  (2)  P  can  limit  rates  of  detritus 
decomposition  and  aufwuchs  production  (Elwood  et  al.  1981)  and  thus  may  control  rates  of 
energy  transformations. 

The  two  sampling  dates  are  representative  of  different  seasons.  On  August  4,  flow  was 
on  the  falling  limb  of  the  summer  storm  hydrograph  (2.0  mVs)  and  coarse  particulate  matter 
(CPM)  transport  was  also  declining  (27.1  mg/m3  at  site  1).  Air  and  water  temperatures  were 
10°C  and  13°C,  respectively.  Flows  were  much  lower  on  October  4  (0.2  mVs)  and  it  had  been 
two  weeks  since  the  last  rain  storm.  Stream  transport  of  autumn  litterfall  (primarily  spruce 
needles)  was  high;  concentration  of  CPM  transported  past  site  1  was  16.5  mg/m3  compared 
with  3.0  mg/m3  of  CPM  transported  past  site  1  on  September  22.  Air  and  water  temperatures 
were  3°C  and  4°C,  respectively. 

On  both  dates,  there  were  no  significant  changes  in  FPP  or  TDP  concentrations  as 
streamwater  was  transported  downstream  (Table  5.1).  However,  both  FPP  and  TDP 
concentrations  were  lower  in  October  than  August  (n=  12,  f=14.1,  P<0.05 ;  n=12,  t=13.0. 
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P<0.05,  respectively)  probably  due  to  lower  flows  and  lack  of  suspension  of  fine  particulate 
materials.  As  well,  in  August,  there  was  no  significant  difference  between  FPP  and  TDP 
concentrations  (n=12,  t= 2.0,  />>0.05)  but  in  October,  TDP  was  significantly  higher  than 
FPP  concentrations  (n= 12,  t=10.0,  F<0.05). 

The  variability  in  P  concentration  between  seeps  was  a  surprising  pattern  that  emerged, 
perhaps  reflecting  different  origins  of  the  seeps  (e.g.  groundwater  vs.  bog)  and  distance 
travelled  before  entering  Two  Creek.  An  overall  trend,  however,  is  relatively  low  FPP 
concentration  and  relatively  high  TDP  concentration  in  seepwater.  The  low  FPP  concentration 
in  seeps  had  no  apparent  effect  (i.e.  no  dilution  effect)  on  stream  water  due  to  the  low  volume 
of  water  entering  the  streams  via  seeps  or  due  to  a  conversion  of  dissolved  to  fine  particulate 
material  as  reported  by  Meyer  and  Likens  (1979). 

These  results  suggest  seasonal  but  no  longitudinal  trends  in  Two  Creek  based  on 
measurement  of  materials  in  transport  due  to:  (1)  significant  differences  in  both  TDP  and 
FPP  concentration  between  August  and  October,  and  (2)  no  downstream  change  in  either  TDP 
or  FPP  concentration  in  August  or  in  October.  Thus,  I  predict  no  change  in  frequency  and 
abundance  of  functional  feeding  groups  of  invertebrates  along  Two  Creek.  However,  I  predict  a 
seasonal  shift  from  collector -filterers  in  August  (Cummins  1974)  to  a  community  comprised  of 
bacteria  and  other  decomposers  in  October.  These  predictions  ignore  the  availability  of  benthic 
food.  However,  these  results  suggest  that  invertebrates  do  not  have  a  large  role  in  transport 
dynamics. 

The  first  predictions  that  functional  groups  would  not  change  in  a  longitudinal  gradient 
appear  congruent  with  trends  observed  in  low -order  streams  in  the  southern  Appalachian  region 
where  macroinvertebrates  accelerated  the  turnover  of  particulate  matter  so  that  by  the  end  of 
the  summer  there  was  almost  no  particulate  material  left  in  the  stream  (Webster  1983). 
However,  I  postulate  that  in  Two  Creek  the  greater  amount  of  particulate  material  transported 
during  summer  as  compared  with  autumn  was  due  primarily  to  higher  flows  and  not 
invertebrate  processing.  I  suggest  that  in  Two  Creek,  P  transport  and  processing  were  more 


»<  K 


* 


56 


similar  to  patterns  observed  at  Hubbard  Brook  (New  Hampshire)  where  macroinvertebrates  do 
not  appear  to  have  a  significant  role,  and  large  amounts  of  detritus  accumulated  and  were 
susceptible  to  transport  during  major  storms  (Fisher  and  Likens  1973,  Meyer  and  Likens  1979). 
The  basis  for  this  suggestion  is  the  strong  pulsing  of  fine  particulate  material  during  storm 
events  in  Two  Creek  (Chapter  2)  which  is  different  from  Hurricane  Branch  where  strong 
pulsing  of  particulate  matter  was  absent  (Webster  1983). 

Thus,  I  suggest  that  abiotic  rather  than  biotic  factors  have  a  dominant  influence  on  the 
partitioning  of  P  transport.  To  either  support  or  refute  these  predictions,  the  following  work 
needs  to  be  done:  (1)  quantifying  the  volume  as  well  as  the  P  concentration  of  water  entering 
the  main  stream  channel  via  seeps,  (2)  quantifying  abundance  of  invertebrates  and  categorizing 
them  according  to  functional  feeding  group,  (3)  measuring  standing  stock  of  detrital  material 
during  all  seasons,  (4)  quantifying  rates  of  invertebrate  detrital  processing  during  all  seasons, 
and  (5)  addition  and  removal  experiments  of  both  invertebrates  and  detrital  material. 
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Table  5.1.  Fine  particulate  phosphorus  (FPP)  and  total  dissolved  phosphorus  (TDP) 
concentrations  (mg/m3)  of  streamwater  at  twelve  different  sites  along  Two  Creek:  site  1  is 
the  main  study  site. 


AUGUST  4  OCTOBER  4 


SITE 

FPP 

TDP 

FPP 

TDP 

1 

10.9 

12.6 

3.9 

9.8 

2 

13.9 

11.7 

3.9 

10.1 

3 

seep 

2.5 

19.5 

7.8 

10.1 

4 

12.4 

11.9 

4.7 

9.4 

5 

11.4 

12.6 

6.6 

8.8 

6 

11.3 

12.2 

5.9 

8.9 

7 

seep 

1.8 

8.7 

0.6 

7.8 

8 

10.9 

12.3 

5.7 

8.3 

9 

11.5 

12.2 

5.8 

8.0 

10 

seep 

4.5 

10.0 

3.2 

6.7 

11 

11.3 

12.3 

5.2 

8.2 

12 

10.8 

12.6 

4.3 

8.9 
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Fig.  5.1.  Two  Creek  study  site  with  numbers  indicating  sampling  sites. 


VI.  APPENDIX  C:  Figures  of  Sakwatamau  River  and  Two  Creek  data  not  included  in  Chapter  II 
and  cluster  analysis  dendrogram  for  total  phosphorus  loading  in  Two  Creek. 
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PHOSPHORUS  CONCENTRATION  (mg/m3) 
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Fig.  6.1.  Discharge,  fine  particulate  phosphorus  (FPP).  and  total  dissolved  phosphorus 
(TDP)  concentrations  at  2-m  intervals  across  the  width  of  Sakwatamau  River. 
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Tooi  ^'r“  ya[ialion  'n  Pr?c'P'tal'on  (PPT),  discharge  (Q),  total  phosphorus  (TP),  and  total  dissolved  phosphorus  (TDP)  concentrations  from  April  to  October, 


PHOSPHORUS  CONCENTRATION  (mg/m3) 
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AUG  5 

AUG  12 

AUG  30 

SEPT  6 

SEPT  14 

SEPT  20 

OCT  14 

OCT  18 


Fig.  6.3 


Fractionation  of  the  fine  particulate  phosphorus  size  fraction  for  Two  Creek. 
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Fig.  6.4.  Fractionation  of  the  fine  particulate  phosphorus  size  fraction  for  Sakwatamau 
River. 
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Fig.  6.5.  Relationship  between  fine  particulate  phosphorus  (FPP)  concentration  and  discharge 
(Q)  for  Sakwatamau  River  (n  =  170).  R  is  rate  of  change  of  discharge. 
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Fig.  6.6.  Relationship  between  total  dissolved  phosphorus  (TDP)  concentration  and  discharge 
(Q)  for  Sakwatamau  River. 
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Fig.  6.7.  Relationship  between  coarse  particulate  matter  (CPM)  and  discharge  (Q)  for 
Sakwatamau  River. 
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Fig.  6.8.  Cluster  analysis  dendrogram  indicating  degree  of  similarity  in  total  phosphorus 
loading  between  months  for  Two  Creek. 


VII.  APPENDIX  D:  Physical  and  chemical  data  collected  from  Two  Creek  and  Sakwatamau 

River  during  1983 
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Table  7.1.  Precipitation  data  collected  from  April -September,  1983,  at  Eagle  Lookout  Station 
(Alberta  Forest  Service  firetower)  located  on  the  northern  edge  of  Two  Creek  watershed. 


Month 

Daily  Maximum 
(mm) 

Average  Daily 
(mm/day) 

Total 

(mm) 

April 

5.0 

0.25 

7.6 

May 

14.7 

1.58 

49.1 

June 

67.3 

8.55 

256.1 

July 

38.9 

6.84 

212.1 

August 

25.4 

1.23 

38.1 

September 

11.0 

1.98 

59.5 

Total 

622.9 

Table  7.2.  Summary  of  1983  discharge  data  for  Two  Creek  and  Sakwatamau  River;  X  is  average  discharge  (mVs)  for  the  sampling  period  and  T  is  total 
discharge  (dam3)  for  the  sampling  period. 
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Table  7.3.  Two  Creek  and  Sakwatamau  River  pH  data  collected  during  1983. 


Date 

Sakwatamau  River 

Two  Creek 

May  14 

8.0 

8.0 

May  27 

8.3 

8.3 

July  31 

7.6 

7.6 

September  6 

8.4 

8.4 

September  22 

8.0 

8.2 

October  1 

8.2 

8.1 

October  14 

8.3 

Mean  (SD) 

8.1  (.29) 

8.1  (.26) 

*  < 


73 


Table  7.4.  Average  nitrate  and  nitrite  nitrogen  (mg/m3),  average  Kjeldahl  nitrogen  (mg/m3), 
and  average  total  nitrogen  (TN;  mg/m3  concentrations  as  well  as  total  nitrogen  to  total 
phosphorus  ratios  (TN:TP)  for  Sakwatamau  River  and  Two  Creek  during  1983. 


Date 

Stream 

Nitrate  +  Nitrite 

Kjeldahl 

TN 

TN :  TP 

May  27 

Sakwatamau  River 

1.9 

458.2 

460.1 

24 

Two  Creek 

1.3 

495.3 

496.6 

27 

August  1 

Sakwatamau  River 

2.4 

603.4 

605.8 

22 

Two  Creek 

3.2 

766.1 

769.3 

26 
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Table  7.5.  Cation  (sodium,  potassium,  calcium,  and  magnesium)  and  anion  (total  alkalinity  as 
CaC03,  bicarbonate,  chloride,  and  sulfate)  concentrations  in  mg/L  from  samples  collected 
from  Sakwatamau  River  and  Two  Creek  on  August  3,  1983. 


Ion 

Sakwatamau  River 

Two  Creek 

Sodium 

6.4 

5.7 

Potassium 

0.2 

0.2 

Calcium 

34.1 

32.1 

Magnesium 

8.9 

8.1 

CaC03 

74.3 

65.6 

Bicarbonate 

90.6 

80.0 

Chloride 

0.10 

0.09 

Sulfate 

0.7 

0.6 

•. »  T : 
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Table  7.6.  Maximum  instantaneous  discharge  (Q;  m3/s),  mean  daily  discharge  (m3/s),  and 
total  annual  discharge  (dam3)  for  an  eleven  year  period  (1973-1983)  collected  by  Environment 
Canada  at  the  mouth  of  Sakwatamau  River  near  Whitecourt,  Alberta. 


Year 

Maximum  Q 

Mean  Daily  Q 

Total  Q 

1973 

65.4 

June  17 

8.1 

174000 

1974 

80.1 

April  30 

9.6 

205000 

1975 

101.0 

June  28 

6.8 

144000 

1976 

58.3 

August  18 

9.6 

204000 

1977 

101.0 

May  30 

1.4 

231000 

1978 

87.8 

June  15 

7.2 

155000 

1979 

14.2 

301000 

1980 

34.5 

May  28 

7.1 

151000 

1981 

50.1 

June  1 

4.7 

98600 

1982 

51.6 

August  4 

6.8 

143000 

1983 

92.6 

June  26 

7.4 

164000 

Mean 

7.6 

179145 

Table  7.7.  Sakwatamau  River  physical  and  chemical  data.  Twenty-one  variables  were  measured  and  are  defined  as  follows: 
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Table  7.8.  Two  Creek  physical  and  chemical  data.  Twenty-one  variables  were  measured  and  are  defined  as  follows: 
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